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ABSTRACT 


!  1  1 


Energetic  chemical  6pecies  have  been  investigated  by  both  theoretical 
and  experimental  methods.  The  spatial  distribution  of  trapped  radicals  in 
y-irradiated  ethyleneglycol  dimethacrylate  polymers  and  of  trapped  electrons 
in  Y*irradiated  methanol,  ethanol  and  2-methyltetrahydrofuran  has  been 
studied  by  paramagnetic  relaxation  methods.  The  radius  within  which  elec** 
trons  are  trapped  increases  with  decreasing  polarity  of  the  matrix.  Charge 
storage  in  y-irradiated  polymethylmethacrylate  (PMMA)  and  other  organic 
polymers  doped  with  electron  traps  has  been  studied  by  measuring  thermally- 
induced  thermoelectric  and  ohmic  currents.  The  conduction  mechanism  for 
the  thermoelectric  current  in  PMMA  appears  to  be  hole  hopping  between 
rotating  carbomethoxy  side  groups.  The  origin  of  the  ohmic  currents  is  not 
yet  clear.  The  theoretical  technique  for  obtaining  lower  energy  bounds  has 
been  extended  to  multi -electron,  polyatomic  molecules  by  using  Gaussian 
basis  orbitals.  The  use  of  Gaussians  allows  all  integrals  to  be  evaluated  in 
closed  form  or  by  a  single  numerical  integration.  In  the  geminal  theory  of 
chemical  bonding  a  new  method  for  the  determination  of  the  best  separated 
pair  geminals  of  a  closed  shell  molecule  has  been  developed,.  The  efficiency 
of  triplet-triplot  energy  transfer  in  naphthalene-biacetyl  mixtures  has  been 
studied. 
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STRUCTURE  AND  REACTIVITY  OF  ENERGETIC  CHEMICAL  SPECIES 


SCOPE  OF  THE  PROGRAM 

The  physics  and  chemistry  of  highly  energetic  chemical  species  such  as 
ions,  electronically  and  vibrationally  excited  molecules,  and  free  radicals 
is  still  in  its  infancy.  The  greatest  paucity  of  information  concerns 
excited  state  species;  their  structure,  reactivity,  and  metastability 
represent  a  newly  developing,  experimentally  difficult,  but  potentially 
rewarding  field  of  study.  Theoretical  techniques  are  just  now  arriving  at 
the  point  at  which  they  can  provide  accurate  and  detailed  information  on 
energetic  species  of  chemical  interest.  Thus,  theory  is  reaching  the  point 
where  it  can  supplement  and  guide  the  experimental  investigation  of 
energetic  species.  The  characteristics  of  energetic  species  of  all  types 
in  the  solid  state  is  also  an  area  of  great  importance  and  limited 
information. 

The  areas  of  study  are  vectored  toward  an  enhanced  understanding  of 
the  physical  and  chemical  characteristics  of  highly  energetic  chemical 
species.  A  balanced  experimental  and  theoretical  study  of  both  structure 
and  reactivity  of  carefully  selected  species  and  systems  utilizing  a 
variety  of  techniques  characterizes  the  general  approach  used.  In  this 
arnual  report  progress  is  reported  in  the  following  areas:  (a)  the  spatial 
distribution  of  trapped  radicals  and  electrons  produced  by  y-irradiation, 
(b)  r ad i at  ion- induced  charge  storage  in  organic  polymers,  (c)  lower  energy 
bound  calculations  on  polyatomic  systems,  (d)  geminal  theory  of  chemical 
bonding,  and  ie)  triplet -triplet  energy  transfer. 


I .  RAD I  AT  I  ON- 1 NDUCED  ENERGY  STORAGE 
IN  ORGANIC  SOLIDS 
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High  energy  radiation  such  as  gamma  irradiation,  vacuum  ultraviolet 
photons  and  solar  wind  protons  expends  its  energy  in  an  organic  solid  by 
ionization  and  excitation  processes.  The  product  ions,  free  radicals  and 
excited  molecules  are  inherently  energetic  chemical  species  but  in  suitable 
solid  matrices  some  of  these  species  can  be  trapped  and  stored.  The  present 
program  has  been  vectored  toward  the  storage  potential,  spatial  distribution, 
and  detrapping  mechanisms  of  trapped  ions  and  free  radicals.  Oir  primary 
focus  has  been  on  polymeric  systems  since  the  trapped  energetic  species  can 
be  "stored"  at  room  temperature  and  above  and  on  ions,  including  electrons, 
since  ions  are  potentially  more  energetic  than  free  radicals. 

A.  SPATIAL  DISTRIBUTION  OF  TRAPPED  RADICALS 
IN  y-  IRRADIATED  ETHYLENEGLYCOL  D I  METHACRYLATE  POLYMERS 

1.  INTRODUCTION 

Radiation  energy  is  deposited  by  y-rays  in  condensed  systems 
inhomogeneously  and  produces  an  initial  inhomogeneous  distribution  of 
energetic  ionized  and  excited  species.  For  very  fast  chemical  reactions 
the  initial  spatial  inhomogeneity  of  the  species  involved  must  be  taken 
into  account.  In  liquid  water  radiolysis  such  reasoning  has  led  to  the 
development  of  the  "spur"  model  and  the  consequent  successful  application 
of  diffusion  kinetics.  Likewise  in  condensed  systems  such  as  polymers 
and  frozen  solutions  it  is  of  considerable  interest  to  determine  if 


3 

trapped  species  exhibit  an  inhomogeneous  spatial  distribution  because  this 
may  constitute  the  factor  that  limits  the  maximum  concentration  of  trapped 
species. 

Recently  we  showed  how  paramagnetic  relaxation  experiments  can  give 

new  insights  into  the  spatial  distributions  of  trapped  electrons  and 

2-5 

hydrogen  atoms  in  frozen  aqueous  systems.  Trapped  electrons  in  irradiated 

alkaline  ices  at  77°K  show  no  change  in  their  relaxation  time  over  a  range 

of  y~dose  from  0.2  to  3  Mrad  even  though  the  electron  concentration  increases 

linearly  in  this  range;  this  implies  that  the  average  spin-spin  interaction 

and  the  local  concentration  of  electrons  remains  constant  in  this  dose 

range.  The  results  are  explained  by  an  inhomogeneous  distribution  of 

electrons  trapped  in  spurs  in  which  only  intraspur  spin-spin  interactions 

are  important.  *  At  higher  doses  the  spurs  overlap,  interspur  spin-spin 

interactions  become  important  and  the  relaxation  time  decreases  with  dose. 

In  contrast,  trapped  hydrogen  atoms  in  irradiated  acidic  ices  at  77°K 

exhibit  a  decreasing  relaxation  time  and  a  linearly  increasing  concentration 

4,5 

over  a  radiation  dose  range  of  0.3  to  5  Mrad.  This  is  the  behavior 
expected  for  a  homogeneous  distribution.  The  trapped  hydrogen  atoms  are 
trapped  only  near  oxyanion  molecules  and  are  expected  to  have  a  uniform 
spatial  distribution.  The  relaxation  results  confirm  this. 

In  the  present  work  we  have  extended  our  paramagnetic  relaxation  method 
to  trapped  organic  free  radicals.  This  method  depends  on  microwave  power 
saturation  and  can  be  applied  only  to  organic  radicals  that  saturate  within 
the  range  of  available  microwave  power  (typically  200  mw).  We  have  found 
that  radicals  trapped  in  irradiated  ethyleneglycol  dimethacrylate  (EDMA) 
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gels  at  room  temperature  fulfill  this  requirement.  The  E£MA  monomers  have 
the  following  structure 

CH3  ch3 

CH2=C“C—0  ( — CH2”CH2“0_i  C-^CHj 
II  V  /n  ll 

o'  0 

and  n  =  1,2,3  and  4  have  been  studied.  Irradiation  produces  a  rigid 

crosslinked  polymer.  Trapped  radicals  in  irradiated  polymethylmethacrylate 

which  is  not  crosslinked  are  also  satur?ible  and  have  been  studied  recently 

6  7 

with  similar  aims  by  Bullock  et  al.  *  The  identity  of  the  EDMA  radicals 
and  the  effect  of  radiation  dose  and  temperature  on  their  relaxation  time 
are  reported.  In  addition,  the  local  radical  concentration  and  the  sample- 
avorage  radical  concentration  have  been  measured  and  compared.  The  results 
indicate  that  the  radicals  in  EDMA  gels  are  of  the  methacrylate  type  and 
are  trapped  with  an  inhomogeneous  spatial  distribution  within  radiation- 
produced  spurs.  The  experimental  results  were  reported  in  more  detail  last 
year  in  Report  No.  AFRPL-TR-67-321 ;  but  the  analysis  was  only  completed 
this  year. 

2.  THEORY 

g 

Port  is  has  classified  the  types  of  paramagnetic  line  broadening  as 
homogeneous  or  inhomogeneous.  An  EPR  line  is  homogeneously  broadened  by 
interactions  which  allow  the  spin  system  to  remain  in  thermal  equilibrium 
during  resonance  absorption.  Inhomogeneous  line  broadening  interactions 
with  nuclear  spins  cause  different  sets  of  electron  spins  to  see  different 
net  local  magnetic  fields.  Each  such  set  of  electron  spins  forms  a  "spin 
packet"  and  the  observed  spectrum  of  the  total  spin  system  is  an  envelope 


MM 
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of  the  superimposed  spin  packets*  The  shape  of  the  envelope  depends  on 
the  intensity  distribution  of  local  magnetic  fields  seen  by  the  electron 
spins  and  the  ratio  of  spin  packet  width  to  the  envelope  width.  If  the 
spin  packet  width  is  much  less  than  the  envelope  width,  as  is  often 
observed,  and  if  the  spin  packets  do  not  interact  magnetically,  then  the 
envelope  shape  is  dependent  only  on  the  intensity  distribution  of  local 
magnetic  fields  and  is  Gaussian. 

An  electron  spin  system  saturates  when  the  population  ratio  between 
the  upper  and  lower  spin  energy  states  deviates  from  its  thermal  equilibrium 
value.  Power  saturation  is  most  readily  studied  by  obtaining  saturation 
curves  vrfiich  are  plots  of  EPR  signal  intensity  versus  microwave  magnetic 
field,  H|,  A  study  of  saturation  behavior  allows  one  to  determine  the 
product  TjT2  or  sometimes  T,  alone;  Tt  is  the  spin-lattice  relaxation  time 
and  Ta  is  the  spin-spin  relaxation  time.  The  relevant  theory  for  obtaining 
relaxation  times  from  saturation  curves  under  slow  passage  conditions  has 
been  summarized  in  a  previous  paper. ^ 

A  typical  saturation  curve  for  the  central  EPR  line  of  free  radicals 

in  tri ethyleneglycol  di methacrylate  gel  is  shown  in  Figure  1-1.  An  ideal 

inhomogeneously  broadened  line  would  have  resulted  in  a  saturation  curve 

which  had  a  slope  of  0  at  high  values  of  H1  whereas  an  ideal  homogeneous 

line  would  have  yielded  a  curve  which  bends  down  at  high  H,  values  even 

3 

more  sharply  than  the  free  radical  curve.  If  homogeneous  saturation 
behavior  is  assumed  T,T2  can  be  calculated  from  Equation  (1) 


V*  =  VrW/2 


(D 


where  y  *®  *h®  gyromagnetic  ratio  of  the  electron  and  =  1.76  x  107  gauss  1  sec 
and  Ht/2  *•  the  value  of  H,  at  which  the  EPR  signal  intensity  is  one-half  of 
what  it  would  have  been  in  the  absence  of  saturation.  If  ideal  inhomogeneous 
saturation  behavior  is  assumed,  TtT2  is  calculated  by  Equation  12). 


T,T*  =  3/y2H?/2  (2) 

For  an  intermediate  case  such  as  that  showi  in  Figure  1-1,  Castner's 
q 

analysis  may  be  applied  using  curve  fitting  procedures  to  correct  H?^2  and 
then  using  Eqn.  (1)  to  obtain  T,T2.  This  method  also  allows  the  calculation 
of  T2  separately  from  Equation  (3) 


wh«re  AtP  is  the  measured  linewidth  at  maximum  slope  of  the  observed 
ms 

Gaussian  line  and  a  is  a  measure  of  the  ratio  of  the  Lorentzian  spin  packet 

width  to  the  observed  Gaussian  linewidth;  a  is  determined  from  the  satura- 

..  9 

tion  curve. 


J 

| 


3.  EXPERIMENTAL 

Mono-,  di-,  tri-  and  tetraethyleneglycol  dimethacrylate  monomers 
(abbreviated  mono—,  di— ,  tri-  ano  tetra^DMAX  where  the  X  indicates  that 
inhibitor  has  been  removed)  were  obtained  from  Sartomer  Resin  Company  and 
contained  60  ppm  hydroquinone  as  an  inhibitor  of  thermal  polymerization. 
The  hydroquinone  was  removed  by  forming  its  sodium  salt  in  1.5  M  NaOH  and 
extracting  with  water.  The  monomer  was  then  dried  over  MgS04  and  stored 
in  a  refrigerator. 
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Samples  were  prepared  by  pipetting  about  2  cc  of  each  monomer  into 
3  mm  o.d.  x  5  cm  Spectrosil  quartz  tubes.  These  were  irradiated  at  about 
35°C  in  a  60Co  Y“*rradiator  at  a  dose  rate  of  0.54  Mrad/hr  as  determined 
by  ferrous  sulfate  dosimetry.  Polymerization  occurred^* ^  and  radicals 
were  trapped  in  the  monomer  gels  which  appeared  to  be  stable  for  at  least 
several  days  after  irradiation. 

The  EPR  measurements  were  made  at  room  temperature  by  placing  the 
irradiated  samples  in  a  quartz  EPR  dewar.  After  a  sample  had  been  measured 
at  room  temperature  it  could  be  measured  at  77°K  by  simply  fill ing  the 
dewar  with  liquid  nitrogen.  All  EPR  measurements  were  made  on  a  Varian-4500 
spectrometer  equipped  with  a  V-4500-41A  L>-Hi  power  microwave  bridge,  an 
audio  modulation  system,  and  a  Varian  dual  cavity  with  dewar  inserts^  The 
bridge  was  operated  in  the  low-power  mode  which  allowed  a  microwave  power 
variation  over  a  40-d8  range.  Slow  passage  conditions  were  achieved  by 

operating  with  a  field  modulation  frequency  of  40  cps  and  a  typical 

« 

modulation  amplitude  cf  0.3  gauss.  The  slow  passage  progressive  saturation 

3 

measurements  were  made  as  described  previously. 

4.  RESULTS 

A  typical  EPR  spectrum  of  radicals  trapped  in  irradiated  tri-£HlAX  is 
shown  in  Figure  1-2;.  This  is  a  direct  recorder  trace  from  the  present  work 
and  shows  the  first-derivative  absorption  signal  obtained  with  fast-passage 
(100  kc)  operating  conditions.  Spectra  from  mono-,  di-,  and  tetra-EBMAX 
are  not  shown  because  they  are  identical  to  that  of  Figure  1-2. 

The  relaxation  time  (TjTj)1/2,  of  free  radicals  trapped  in  tri-£EMAX 
was  measured  as  a  function  of  radiation  dose  and  temperature  under 


slow-passage  conditions.  The  central  line  of  the  EPR  spectrum  was  used  for 

the  measurements.  From  saturation  curves  such  as  the  one  shown  in  Figure  1-1 

the  value  of  was  measured  and  (T^)1/2  was  calculated  from  Equation  (2) 

q 

assuming  ideal  inhomogeneous  conditions.  Castner's  method  was  also  used 
to  calculate  (T^)1/2  [Equation  0)3  and  T2  [Equation  (3)].  These  data  ai*e 
•presented  in  Table  l-l . 

Table  l-l I  contains  measurements  of  the  linewidth  between  derivative 
maxima,  Ah^s,  of  the  central  EPR  line  versus  for  trapped  free  radicals 
in  tri-EDMAX  at  room  temperature. 


5.  DISCUSSION 

a.  B&dical  .Identificalign 

The  9-line  EPR  spectrum  of  Figure  1-2  for  radicals  trapped  in  tri-EDMAX 

is  similar  to  those  for  radicals  in  mono-,  di-,  and  tetra-EDMAX.  The  only 

observable  difference  among  them  is  that  the  trapped  radical  yield  increases 

in  the  direction  mono*di+tr i-*  tetra-EDMAX.  Relative  yield  ratios  for  the 

mono-  through  tetra-EDMAX  monomers  are  1.0:1. 3:2. 5*2. 6.  The  increase  in 

radical  yield  correlates  with  the  increase  in  length  of  monomer  molecule 

and  density  of  monomer  solution.  It  is  concluded  that  the  radical  species 

is  the  same  in  all  four  monomer  gels. 

Irradiated  polymethylmethacrylate  IPMMA)  produces  a  radical  with  an 

EPR  spectrum  identical  to  that  in  the  EDMAX  gels.^  This  spectrum  has 

received  a  great  deal  of  study  and  there  has  been  considerable  argument 

13 

as  to  whether  the  9-line  spectrum  is  due  to  one  or  two  radicals.  However, 

the  model  suggested  by  Symons^  and  supported  by  the  high  resolution 

15 

experiments  of  Fischer  and  the  deuteration  experiments  of  Kourim  and 
16 

Vacek  lead  us  to  assign  the  spectrum  in  Figure  1-2  to  the  single  radical 
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Table  I -I  I 

EPR  Linewidth  of  Free  Radicals 
in  Tri-EDMAX  Versus  Microwave  Magnetic  Held* 


Gauss 

Hi  Gauss 

5.2 

0.00354 

5.2 

0.0121 

5.2 

0.0238 

5.2 

0.0900 

Measured  on  Central 

Line  of  EPR  Spect 
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CH, 

r-ch2-c* 

I 

COO  (CHjCHjO)  3C0C  (CH3)=CHj 

CHj 

where  R  is  the  polymer  chain.  The  quartet  from  the  CH3  group  (a  =  23  G) 
it  split  by  the  two  nearly  equivalent  methylene  protons  (a  2  *=  11  and  14  G)  3 
to  give  the  observed  9-line  spectrum;  resolution  of  the  even-numbered  spectral 
lines  is  not  achieved. 

b.  Line  Shape  aq4  Broaden img  MgchflOilB 

Table  l-ll  shows  that  the  central  line  of  the  monomer  radical  EPR 

spectrum  does  not  widen  appreciably  at  the  highest  microwave  powers  available 

from  the  klystron.  The  average  line  shape  factor  for  this  line  was  measured 

to  be  2.3;  2.2  is  the  theoretical  value  for  a  pure  Gaussian  line.  These 

two  results  imply  that  the  EPR  line  is  broadened  mainly  by  inhomogeneous 

mechanisms.  The  broadening  is  principally  due  to  unresolved  hyperfine 
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interaction  and  to  the  nuclear  spin  state  of  the  p-protons.  Analysis 
of  progressive  saturation  curves  (Figure  1-1 )  by  Castner’s  method  results 
in  a  value  of  approximately  0.5  for  his  "a"  parameter  [Equation  (3)].  Although 
Castner's  method  involves  considerable  uncertainty  when  "a"  is  as  large  as 
0.5»  it  does  indicate  that  the  spin  packet  width  is  not  negligible  in 
comparison  to  the  total  linewidth. 

C.  Dose  aq4  Temperature  Effects  £Q  Spatial  Pi atfifaut LPDa  SiL  Radicals 

Table  l-l  summarizes  relaxation  time  as  a  function  of  radiation  dose 
and  of  temperature  for  radicals  trapped  in  tri-EDMAX  gels, 
been  determined  by  Portis'  method  as  well  as  by  Castner's  method.  Castner's 


a 
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method  has  also  bean  applied  to  the  EPR  central  linewidths  to  calculate  Ta 
(Equation  3)  separately.  The  data  show  that:  (a)  (TjT*)1/*  and  Ta  do  not 
change  significantly  with  increasing  radiation  dose  to  2  Mradj  and 
(b)  Ta  appears  to  increase  upon  going  from  29S°K  to  77°K,  whereas  Tt  does 
not  change  significantly. 

Bullock  and  Sutcliffe^  have  made  similar  measurements  on  methyl¬ 
methacrylate  radicals  and  it  is  interesting  to  compare  the  results  obtained 
in  this  work  for  tri-EDMAX  radicals  which  have  a  similar  structure  with 
their  results  for  radicals  trapped  in  photopolymer i zed  pure  methylmethacrylate. 
They  obtain  a  value  of  2.9  x  10*6  sec  at  299°K  for  (TiT*)1/*  by  Castner's 
method  which  compares  favorably  with  the  value  3.1  x  10”6  obtained  in  the 
present  work.  Their  values  of  5.7  x  10"*  sec  for  Ta  and  2.5  x  10  *  sec  for 
Tt  are  in  fair  agreement  with  Ta  ~  3*6  x  10"*  sec  and  T1  =  2.7  x  10  *  sec 
obtained  here.  Their  value  of  Castner's  tta"  parameter  Z  0.30  and  considering 
the  precision  of  the  method  is  comparable  to  the  value  "a"  Z  0.5  measured 
in  this  work.  It  is  clear  from  these  comparisons  that  the  relaxation 
characteristics  of  trapped  methylmethacrylate  radicals  are  very  similar 
to  those  for  the  trapped  radicals  in  tri-EEMAX  gels.  The  EDMAX  gels  are 
crosslinked  while  PMMA  is  notf  and  the  mobility  of  segments  is  lower  in  the 
EDMAX  gels.  These  structural  features  apparently  do  not  affect  the 
relaxation  characteristics. 

From  0.5  to  2.0  Mrad  Ta  remains  constant  although  the  concentration 

of  radicals  is  increasing  linearly  in  this  dose  range.  In  analogy  to  our 
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previous  interpretations  this  indicates  that  the  radicals  are  trapped 
in  spurs  and  that  the  local  radical  concentration  is  higher  than  the  average 
radical  concentration.  These  organic  radicals  exhibit  the  same  relaxation 


15 


behavior  as  do  trapped  electrons  in  alkaline  ices.  At  higher  doses  the 
valuo  of  Ta  should  decrease  due  to  overlap  ;o’f  spurs. 

The  value  of  Castner's  "a"  allows  one  tb  calculate  the  dipolar  spin 
packet  linewidth*  from  Equation  (4). 


ift  19 

Vyard  has  modified  Kittel  and  Abrahams'  treatment  of  dipolar  linewidth 

to  apply  to  non-equivalent  magnetic  centers  in  glassy  and  polycrystalline 

samples*  and  obtains  =  32  M  where  M  is  a  local  molarity  of  trapped 

spins  contributing  to  the  dipolar  linewidth.  The  radicals  in  EDMAX  gels 

may  be  equivalent  in  which  case  =  48  11}  we  shall  use  the  value  of 

48  M  to  calculate  the  local  spin  molarity.  The  calculated  local  spin 

molarity  is  0.037  M  and  is  constant  from  0.5  to  2  Mrad.  Assuming 

Q(radical)  -  4  the  average  spin  concentration  at  2  Mrad  is  only  0.0064  M 

which  is  a  factor  of  5  less  than  the  local  spin  molarity.  A  similar 

analysis  on  trapped  radicals  in  polymethylmethacrylate  also  shows  a  ratio 

7 

of  local  to  average  radical  concentration  greater  than  one.  The  actual 
ratio  of  local  to  averse  spin  molarity  at  a  given  dose  depends  on  the 
quantitative  validity  of  -  48  M  and  on  the  precision  of  a*  but  it  is 
certainly  greater  than  one  at  doses  up  to  2  Mrad.  This  is  completely 
consistent  with  the  spur  concept  deduced  from  radiation  dose  effects. 

The  temperature  data  indicate  that  Ta  increases  by  6(#  as  the 
temperature  is  lowered  from  299°K  to  73°K.  This  curious  temperature  de¬ 
pendence  also  is  found  for  radicals  in  irradiated  polymethylmethacrylate/ 
The  increase  in  Ta  indicates  that  the  local  radical  concentration  is 


decreased  at  lower  temperatures.  This  could  be  related  to  contraction  of 
adjacent  polymer  chains*  but  the  magnitude  of  the  effect  is  similar  in  linear 
polymethylmethacrylate  and  in  cross-linked  EDMA  gels.  The  effect  of  a 
wider  temperature  range  on  T2  needs  to  be  studied  before  more  speculation 
on  the  cause  of  the  temperature  dependence  is  warranted. 
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B.  PARAMAGNETIC  RELAXATION  OF  TRAPPED  ELECTRONS 
IN  Y- IRRADIATED  ORGANIC  GLASSES 


1.  INTRODUCTION 

Recently  we  showed  how  paramagnetic  relaxation  experiments  can  give 

useful  and  semi -quaint i tat ive  insights  into  the  spatial  distributions  of 

1-5 

trapped  electrons  and  hydrogen  atoms  in  frozen  aqueous  systems. 

Radiation-produced  electrons  are  also  produced  in  organic  glasses  at  77°K 
6  7-9 

such  as  alcohols  and  2-methyl tet r ahydr of uran  (MTHF)  and  can  be 

detected  by  optical  and  EPR  spectroscopy.  Iffe  report  here  the  paramagnetic 
relaxation  characteristics  of  trapped  electrons  in  these  organic  glasses. 

The  radiation  dose  dependence  of  the  relaxation  time  gives  information 
about  spatial  distributions*  while  the  relaxation  time  magnitudes  and  the 
linewidths  give  information  about  the  nature  of  the  electron  trapping  site. 

A  revealing  correlation  has  been  found  between  the  size  of  the  spatial 
inhomogeneities  associated  with  radiation-produced  trapped  electrons  and  the 
polarity  of  the  matrix. 

2.  EXPERIMENTAL 

Reagent  grade  methanol  *  ethanol  and  MTHF  were  degassed  on  a  vacuum 
line  and  dried  over  a  freshly  formed  sodium  film.  Most  methanol  and 
ethanol  samples  were  prepared  as  -2  mm  glassy  spheres  by  dropping  drops  of 
the  solution  into  liquid  nitrogen;  this  involved  a  short  exposure  to  air. 
MTHF  samples  were 'prepared  by  sealing  into  3  mm  o.d.  Spectrosil  quartz 
ampoules  or  by  making  glassy  spheres  without  exposure  to  air.  The  latter 
was  accomplished  by  pressurizing  the  vacuum  line  with  dry  nitrogen  to 


force  out  drops  through  a  greaseless  stopcock  into  liquid  nitrogen;  the 
short  distance  between  the  stopcock  and  the  liquid  nitrogen  surface  was 
purged  with  dry  nitrogen. 

Irradiations  were  carried  out  at  77°K  in  a  U.  S.  Nuclear  Co60 
y-ir radiator  at  a  dose  rate  of  0.51  Mrad/hr  as  determined  by  ferrous  sulfate 
dosimetry.  For  integral  dose  calculations  the  measured  dose  rate  was 
corrected  for  radioactive  decay. 

The  slow  passage  progressive  saturation  measurements  were  made  as 
described  before  at  73°K.  Bubbling  of  liquid  nitrogen  was  prevented  by 
using  a  helium  gas  stream.  Slow  passage  conditions  were  achieved  by 
operating  with  a  field  modulation  frequency  of  40  hz  and  a  typical  modula¬ 
tion  amplitude  of  0.5  gauss.  Progressive  saturation  curves  were  obtained 
by  normalizing  the  EPR  signal  heights  to  a  constant  amplifier  gain  setting 
and  plotting  these  heights  versus  microwave  magnetic  field,  H,.  H,  was 
measured  as  described  previously.  Each  sample  using  spheres  utilized  3-5 
spheres.  Most  samples  were  run  in  duplicate  or  triplicate. 

3.  RESULTS 

The  EPR  spectrum  of  yirradiated  methanol  at  77°K  consists  primarily  of 
a  triplet.  This  triplet  is  shown  in  Figure  1-3  and  is  due  to  the  CH20H  - 
radical  superimposed  on  a  singlet  due  to  e^~.  Two  small  outer  lines  are 
also  noted  in  the  figure;  these  are  split  by  130  gauss  and  are  due  to  the 
CHO  radical.  Their  yield  is  very  small.  The  CH20H  radical  spectrum  has  a 
theoretical  1;2:1  intensity  ratio.  The  excess  intensity  in  the  center 
line  is  due  to  the  trapped  electron;  this  can  be  demonstrated  by  optical 
bleaching  in  the  visible  of  the  trapped  electron  line.^  The  radical  spec¬ 
trum  was  manually  subtracted  from  the  total  spectrum  to  obtain  the  trapped 


electron  signal  intensity  as  a  function  of  microwave  power.  This  procedure 

led  to  relatively  consistent  results  although  the  accuracy  is  not  high. 

The  trapped  electron  line  in  methanol  had  a  linewidth  at  maximum  slope  if 

11  gauss  and  a  line  shape  that  was  approximately  Gaussian.  Figure  1-4 

shows  a  typical  power  saturation  curve  for  the  trapped  electron  in  methanol. 

The  relaxation  times  were  calculated  from  the  saturation  curves  by  two 

methods.  Method  1  was  the  Port  is  method  in  which  ideal  inhomogeneous 

broadening  behavior  was  assumed.'  "he  second  method  was  Castner's 

analysis  of  the  nonideal  inhomogeneous  broadening  case.^ 

The  EPR  spectrum  of  y irradiated  MTHF  is  shown  in  Figure  1-5.  The 

MTHF  spectrum  consists  of  a  singlet  e^“  line  superimposed  on  a  7-line 

radical  spectrum.  The  7-line  radical  spectrum  has  been  interpreted  as  due 

to  a  free  radical  site  at  the  number  two  carbon  position.  The  hyperfine  is 

then  due  to  three  methyl  protons  and  one  of  the  beta  CH2  protons  with 

splitting  constant  A  and  the  other  beta  CH2  proton  with  splitting  constant 

2A  to  give  relative  intensities  of  1 :4:7:8:7:4:1 •  As  can  be  seen  in 

Figure  1-5  the  e^~  line  is  much  more  intense  than  the  underlying  radical 

spectrum  so  that  the  radical  spectrum  can  be  easily  and  accurately  subtracte: 

from  the  total  spectrum.  Thus  the  results  obtained  for  e^”  in  MTHF  are  much 

more  accurate  than  those  obtained  for  e^”  in  methanol.  The  trapped  electron 

in  MTHF  had  a  linewidth  at  maximum  slope  of  3.8  ±  0.2  gauss  and  had  a  line 
2  12 

shape  parameter  *  of  2.2  which  indicated  that  it  was  closely  Gaussian.  A 

13 

typical  power  saturation  curve  for  e^’  in  MTHF  is  shown  in  Figure  1-6. 

Note  that  the  peak  of  the  saturation  curve  in  MTHF  occurs  at  a  value  of  H! 
five  times  lower  than  in  methanol.  The  MTHF  saturation  curves  were  analyzed 
by  both  Port  is'  and  Castner's  methods. 
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t 


25 


A  few  experiments  were  also  done  on  ethanol.  It  exhibits  a  5-line 
radical  spectrum  superimposed  on  a  singlet  trapped  electron  spectrum. 
Saturation  curves  were  obtained  by  subtracting  the  radical  spectrum  from  the 
total  spectrum  but  the  accuracy  was  not  high.  In  general,  the  results  in 
ethanol  paralleled  those  in  methanol. 

Table  l-lll  summarizes  the  parameters  determined  for  the  trapped 
electrons  in  the  various  matrices.  Note  that  the  relaxation  times  for  the 
electrons  in  the  alcohols  are  about  an  order  of  magnitude  smaller  than  the 
relaxation  times  for  trapped  electrons  in  MTHF  or  in  the  sodium  hydroxide 
ice.  All  of  the  ■results  in  Table  i-l  1 1  are  at  doses  low  enough  to  be  in 
the  dose  independent  region  of  the  relaxation  times. 

Figure  1-7  shows  the  variation  of  relaxation  time  versus  radiation  dose. 
Determinations  of  Ta  show  that  the  dose  variation  in  CT ,T2) ,/a  is  due 
largely  to  variations  in  Ta.  In  Figure  1-7,  however,  values  of  (TjTa)1/2 
determined  by  Port  is'  method  are  plotted  since  they  could  be  determined 
for  a  greater  number  of  samples.  Note  that  the  relaxation  time  for  electrons 
in  methanol  remains  constant  from  0.5  to  2.0  Mrad  and  then  decreases.  In 
contrast,  the  relaxation  time  for  electrons  in  MTHF  remains  constant  only 
from  0.1  to  1.0  Mrad  before  decreasing.  Both  of  these  contrast  with  the 
results  of  trapped  electrons  in  alkaline  ices^  in  which  the  relaxation  time 
remains  constant  to  about  4  Mrad  dose.  Figure  1-8  shows  relaxation  data  for 
electrons  in  methanol  extended  to  very  high  doses  and  also  shows  the  total 
spin  concentration.  The  decrease  in  the  relaxation  time  levels  off  at  the 
same  dose  at  which  the  increase  in  the  total  spin  concentration  reaches  a 
plateau. 
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4.  DISCUSSION 

«.  Una  Shaft*  aad  Brgadsaiog  Mechanism 


The  line  shape  of  e^“  in  liTHF  fits  a  Gaussian  function  quite  well 


according  to  both  our  data  and  the  data  of  others.^  This  line  shape  is 
maintained  at  doses  up  to  6  Mr ads;  at  higher  doses  it  deviates  slightly 
from  a  good  Gaussian  fit.  The  line  shape  of  e^~  in  methanol  and  ethanol  is 
also  Gaussian.  However,  the  fit  cannot  be  determined  precisely  because  of 
overlap  with  the  CH^OH  radical.  The  e^~  EPR  linewidths  in  MTHF,  methanol 
and  ethanol  do  not  increase  appreciably  during  power  saturation.  At  the 
higher  microwave  powers  small  increases  in  the  e^~  linewidths  would  be 
undetectable  because  of  the  increased  size  of  the  radical  signal  relative  to 
that  of  e^~.  However,  it  is  clear  that  the  e^“  lines  do  not  broaden  with 
H,  to  the  extent  theoretically  predicted  for  homogeneously  broadened  lines. 

These  two  results  imply  that  the  e^”  EPR  lines  are  broadened  mainly  by 

inhomogeneous  mechanisms.  The  shape  of  the  power  saturation  curves  in 

* 

Figures  1-4  end  1-6  are  consistent  with  this  conclusion.  Analysis  of  the 
saturation  curves  by  Castner's  method^  give  values  of  his  4  parameter  of 
less  than  one  which  indicate  that  the  spin  packet  width  is  not  negligible 
compared  to  the  total  linewidth.  The  most  likely  inhomogeneous  broadening 
mechanisms  of  importance  are  unresolved  hyperfine  interaction  with  surrounding 
protons  and  anisotropy  broadening.  The  former  mechanism  can  he  evaluated 
semi-quant itatively  from  deuterated  compounds^  The  linewidth  of  e^”  in 
ethanol-dg  is  about  5  gauss;  a  precise  value  cannot  be  obtained  due  to 
interference  by  the  radical  spectrum.  Deuterated  MTHF  was  not  examined  due 
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the  large  linewidth  decrease  in  ethanol-dg  does  indicate  that  unresolved 
hyperfine  interaction  is  probably  the  dominant  broadening  mechanism. 

b.  Dose  Effects  on  Spatial  Distributions  of  e+  and  Radicals 

Studies  in  which  relaxation  time  is  measured  as  a  function  of  radiation 
dose  yield  information  about  interactions  between  unpaired  electron  spins. 

For  a  uniform  distribution  of  trapped  species  the  average  distance  between 
them  decreases  with  increasing  radiation  dose,  and  the  spin-spin  inter¬ 
action  becomes  stronger.  In  this  case  Ta  is  expected  to  decrease  with 
increasing  dose.  This  behavior  is  observed  for  trapped  H  atoms  in  acidic 
ices.^  However,  Figure  1-7  shows  that  T2  for  e^~  in  polar  organic  glasses 
behaves  somewhat  differently.  For  e^~  in  MTHF  and  in  methanol  the  dose 
variation  in  CT,Ta) is  due  to  a  variation  in  Ta.  The  important  result  is 
that  (T  1TJe)  for  e^“  in  methanol  remains  constant  at  5.2  x  10"6  sec  from 
0.5  to  2.0  Mrad  dose  before  decreasing.  Similarly  (T,^)1/2  for  e^”  in 
MTHF  remains  constant  at  3.6  x  10-5  sec  from  0.1  to  1.0  Mrad  before  decreasing. 

The  constancy  of  (T ,Ta) over  a  given  dose  range  indicates  that  e^” 

has  spin-spin  interactions  with  unpaired  spins  very  close  to  it  and  only 

weakly  interacts  with  other  much  more  distant  spins.  This  suggests  that  the 

spins  are  trapped  in  radiation-produced  regions  of  inhomogeneity  which  we 

shall  call  spurs.  At  higher  doses  the  spurs  overlap  and  (T,^)1/2  begins  to 

decrease.  Similar  results  and  interpretation  have  been  reported  for  e^“  in 
2 

alkaline  ice. 

The  inverse  relationship  between  the  var  iation  in  (T1T2)1^2  for  e^“ 
and  the  total  spin  concentration  with  radiation  dose  is  demonstrated  for 
methanol  in  Figure  1-8.  The  total  spin  concentration  includes  e^-  and  CHa0H 


spins;  since  the  EPR  spectra  of  e^  and  CHa0H  overlap  both  types  of  spins 

undergo  relaxation  with  e^  .  Note  that  as  the  total  spin  concentration 

reaches  a  limiting  maximum  value  above  18  Mrad,  the  dose  variation  in  (T,Ta)^ 

also  reaches  a  limiting  minimum  value. 

A  semi --quantitative  measure  of  the  spur  size  may  be  obtained  from  the 

dose  at  vrfiich  (T,Ta)^  begins  to  decrease.  This  dose  is  1  Mrad  in  MTHF, 

2 

2  Mrad  in  methanol  and  about  4  Mrad  in  alkaline  ice.  We  assume  that  at 
these  doses  the  spurs  just  overlap  to  fill  the  entire  volume  and  that  there 
are  G(e^~)  electron  spins  per  spur  where  G  is  the  yield  per  100  ev  of 
absorbed  radiation  energy.  Then  the  spur  radius  in  angstroms  is  given  by 


where  p  is  the  density  and  D  is  the  dose  in  Mrad  at  which  O^Tj)^  begins  to 
decrease.  The  calculated  spur  radii  are  72  A  in  MTHF,  59  ^  in  methanol  and 
42  \  in  alkaline  ice.  This  calculation  leads  to  an  overestimate  of  the  spur 
radii,  but  the  trend  with  matrix  polarity  is  unambiguously  shown  by  the 
experimental  data.  As  the  matrix  polarity  decreases  from  alkaline  ice  to 
methanol  to  MTHF  the  size  of  the  spurs  increases  and  the  amount  of  spatial 
nonuniformity  exhibited  by  the  trapped  radicals  becomes  less.  In  a  nonpolar 
matrix  like  3-methylpentane  glass,  trapped  electrons  produced  by  y-irradiation 
are  expected  to  show  little  if  any  evidence  of  spatial  inhomogeneity.  The 
correlation  between  spur  size  and  matrix  polarity  implies  that  radiation- 
produced  electrons  travel  further  before  being  trapped  the  less  polar  the 
matrix.  This  suggests  that  the  trapping  cross  section  or  the  trap  density  for 
radiation-produced  electrons  is  higher  in  more  polar  matrices. 


v. 

The  dependence  of  (T , T a on  radiation  dose  shows  clearly  but  qualitatively 
that  spatial  inhomogeneity  of  trapped  electrons  and  radicals  exists*  The 
absolute  value  of  T2  or  equivalently}  }  the  dipolar  spin  packet  linewidth} 
allows  one  to  calculate  a  local  spin  concentration  to  compare  with  the 
sample  average  spin  concentration.  The  value  of  T2  and  depend  on  the  a 
parameter  which  is  only  accurate  to  about  3056  for  a  >  0.2*  However,  to  the 
accuracy  of  a»  ^  may  be  calculated  from  Eqn.  (l).^ 


ms 
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ms 
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1*47 


(1) 


The  relation  between  the  dipolar  linewidth  and  the  local  spin  concentration 

was  given  by  Kittel  and  Abrahams^*  for  single  crystals  and  was  modified  by 
15 

Wyard  to  apply  to  non~equi valent  magnetic  centers  in  glassy  and  poly- 
crystalline  matrices*  ifyard  obtains  of  =  3-2  M  where  M  is  the  total 
molarity  of  trapped  spins  contributing  to  the  dipolar  linewidth.  The  local 
molarity  is  then  calculated  from  Eqn.  (2)  and  the  sample  average  molarity 


M( local)  = 


ms 


32  x  1.47 


(2) 


is  calculated  from  Eqn.  (3)  where  0  is  the  dose  in  Mrad  and  p  is  the  density. 


M(sample  average)  =  ?f*0*aLi?Efos)Dp  (3) 


The  results  for  1  Mrad  total  dose  are  summarized  in  Table  i-IV.  For  each 
matrix  M( local )AHaverage)  is  greater  than  one  which  again  indicates  a  non- 
uniform  spatial  distribution  of  trapped  radicals. 


34 


5.  REFERENCES 

1.  J.  Zimbrick  and  L.  Kevan,  J.  Am.  Chem.  Soc.,  8§,»  3678  (1966). 

2.  J.  Zimbrick  and  L.  Kevan,  J.  Chem.  Phys.,  42#  2364  (1967). 

3.  J.  Zimbrick  and  L.  Kevan,  Nature,  214#  693  (1967). 

4.  J.  Zimbrick  and  L.  Kevan,  J.  Chem.  Phys.,  42#  5000  (1967). 

5.  L.  Kevan  and  D.  H.  Chen,  J.  Chem.  Phys.,  42#  1970  (1968). 

6.  C.  Chachaty  and  E,  Hayon,  Nature,  200.  59  (1963). 

7.  M.  R.  Ronayne,  J.  P.  Guurino  and  W.  H.  Hamill,  J.  Am.  Chem.  joc»,  §4, 
4230  (1962). 

8.  G.  A.  Salmon,  Disc.  Faraday  Soc.,  22#  284  (1963). 

9.  D.  Smith  and  J.  J.  Pieroni,  Can.  J.  Chem.  42,  876  (1965). 

10.  A.  M,  Port  is,  Phys.  Rev.,  21#  I0?1  (1953). 

11.  T.  G.  Castner,  Phys.  Rev.,  115.  1506  (1959). 

12.  G.  E.  Pake  and  E.  M.  Purcell,  Phys.  Rev.,  24#  1i®4  (1948). 

13.  Our  curve  is  in  substantial  agreement  with  recent  data  of  J.  Lin, 

K.  Tsuji  and  F.  Williams,  J.  Am.  Chem.  Soc.,  22#  2766  (1968). 

14.  C.  Kittel  and  E.  Abrahams,  Phys.  Rev.,  22#  238  (1953). 

15.  S.  J.  Wyard,  Proc.  Phys.  Soc.,  8£>  587  (1965). 


C.  RADI  AT  I  ON- INDUCED  CHARGE  STORAGE  IN  ORGANIC  POLYMERS 

(with  James  Barnes) 
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1.  INTRODUCTION 

A  distinguishing  effect  of  high  energy  radiation  on  matter  is  that  it 
produces  ionization.  This  can  be  both  desirable  and  deleterious.  If  the 
ions  can  be  trapped  and  stored  they  represent  potential  stored  energy.  The 
efficiency  and  maximum  amount  of  energy  storage  is  partially  dependent  on 
the  spatial  distribution  of  the  trapped  species.  This  aspect  was  investi¬ 
gated  in  previous  sections  A  and  B.  In  this  section  we  concern  ourselves 

*>  i 

with  fundamental  investigations  of  factors  that  affect  the  utilization  of 
stored  energy  in  the  form  of  ions.  These  factors  include  the  detrapping 
mechanism  and  the  conduction  mechanism  of  the  trapped  charges. 

Conduction  mechanisms  of  trapped  charges  in  polymeric  systems*  especially 
pr lymethylmethacrylate  (PMMA),  have  been  investigated  and  represent  a 
continuation  of  the  studies  reported  in  the  1967  Annual  Report.  Thermal 
heating  is  used  to  detrap  the  trapped  charges  in  y-irradiated  polymers 
doped  with  electron  traps.  Conduction  of  these  charges  in  a  preferred 
direction  was  studied  by  application  of  a  temperature  gradient  across  the 
sample  and  by  an  applied  electric  field.  The  characteristics  of  thermo¬ 
electric  currents  generated  by  a  temperature  gradient  in  irradiated  PMMA 
indicate  that  the  positive  charge  becomes  mobile  when  the  carbomethoxy  side 
groups  of  the  polymer  chain  begin  to  rotate  freely.  This  finding  suggests 
that  polymers  with  a  proscribed  3-dimensional  structure  can  be  designed  to 
allow  charge  storage  to  persist  to  quite  high  temperatures.  Currents  generated 


by  applied  electric  fields  are  less  clearly  understood  and  seem  to  be 
independent  of  the  thermoelectric  currents.  Further  work  is  aimed  toward 
clarification  of  this  point. 

2.  EXPERIMENTAL 

Samples  of  commercial  polymers  were  machined  from  cylindrical  rods  into 
discs,  normally  2.5  cm  in  diameter  and  1.0  cm  thick.  Attachment  of  the 
electrodes,  irradiation  procedure,  current  measurement  apparatus  and  other 
experimental  details  have  been  described  in  the  1967  Annual  Report.  Only 
new  experimental  details  will  be  presented  here. 

Synthesized  PMMA,  both  pure  and  with  additives  designed  to  act  as 
electron  or  hole  traps,  was  prepared  by  thermal  polymerization  of  degassed, 
inhibitor-free  methylmethacrylate  (Rohm  and  Haas  Co.).  To  avoid  bubble 
formation  the  synthesized  samples  were  prepared  in  small  tubes  and  the 
typical  sample  size  was  1.0  cm  in  diameter  and  0.4  cm  thick.  The  temperature 
cycle  used  for  the  thermal  polymerization  of  MMA  in  this  work  consisted  of 
an  initial  period  at  90°C,  a  period  of  60°C,  followed  by  heating  at  75°C, 
and  90°C.  The  initial  temperature  was  maintained  until  the  viscosity  of  the 
liquid,  when  hot,  was  like  that  of  glycerine.  The  time  necessary  to  reach 
this  viscosity  varied  with  the  tube  diameter  and  the  length  of  time  the 
monomer  had  been  stored  at  room  temperature  prior  to  use.  This  normally 
ranged  from  1  to  5  hours.  The  6C°C  temperature  cycle  was  maintained  for 
approximately  16  hours  followed  by  6  to  8  hours  at  75°C.  Hardening  at  90°C 
for  ipproximately  16  hours  completed  the  polymerization. 

Figure-  1-9  shows  the  circuit  block  diagram  for  the  measurement  of 
thermally- induced  charge  conduction  in  irradiated  polymers  under  the  influence 
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Figure  1-9.  Circuit  block  diagram  for  current 
measurement  in  an  electric  field. 


of  an  electric  field.  The  irradiated  samples  were  heated  in  the  absence  of 
a  temperature  gradient.  The  electric  field  was  applied  by  connecting  a  DC 
power  supply  between  the  right-hand  electrode  of  the  sample  holder  and  the 
circuit  ground.  The  power  supply  was  a  Keithly  Model  240A  High  Voltage 
Supply  with  an  output  of  0  to  1200  volts.  The  polarity  of  the  electric  field 
will  refer  to  the  polarity  of  the  right-hand  electrode  of  the  sample  holder 
relative  to  ground. 

In  order  to  separate  the  effects  of  charge  flow  between  the  sample 
and  the  external  circuit  from  internal  charge  movement,  measurements  were 

I 

made  with  layers  of  insulating  material  placed  between  the  sample  and  the 
sample  holder  electrodes.  Glass  microscope  slides,  1  mm  thick,  were  cut  into 
squares  (2$  mm  x  25  mm)  and  used  as  "blocking  electrodes".  One  side  of  each 
square  was  made  conducting  by  coating  it  with  Eccobond  solder.  The  tabs  of 
each  sample  to  be  measured  with  blocking  electrodes  were  broken  off  of  the 
surface  of  the  sample  after  irradiation.  The  samples  were  mounted  in  the 
sample  holder  with  a  blocking  electrode  on  either  side.  Each  blocking 
electrode  was  mounted  with  its  conducting  surface  in  contact  with  the  electrode 
of  the  sample  holder.  Current  measurements  were  then  performed. 

The  electron  paramagnetic  resonance  (EPR)  measurements  were  made  with  a 
Varian  V-4500  EPR  system  with  100  kc  magnetic  field  modulation.  Samples  of 
Lucite  and  polyvinyl  chloride  (PVC)  were  cut  from  1/8-inch  diameter  cylind¬ 
rical  rod.  Each  sample  was  approximately  3  inches  long,  and  was  handled  by 
a  thin  thread  tied  to  one  end.  An  absorbed  dose  of  0.5  Mrad  was  given  to 
each  sample  before  it  was  mounted  in  the  EPR  cavity.  During  measurement  the 
samples  were  heated  at  a  linear  rate  of  6°C/min  by  circulating  hot  nitrogen 
gas  through  the  EPR  cavity.  The  gas  was  heated  by  passing  it  over  a  heating 


coil*  The  rate  of  heating  the  sample  was  controlled  by  the  rate  of  gas  flow 
and  the  amount  of  current  flowing  through  the  heating  coil.  The  temperature 
of  the  sample  in  the  EPR  cavity  was  measured  with  a  Digitec  Digital  Thermo- 
couple  Thermometer.  The  thermocouple  was  placed  on  the  surface  of  the  sample 
and  sensed  a  temperature  slightly  higher  than  the  average  temperature  of  the 
entire  sample.  EPR  spectra  were  taken  periodically  throughout  the  heating 
cycle  from  25°C  to  100°C. 

3.  RESULTS 

Most  of  our  results  on  thermoelectric  currents  in  Lucite  including 
dose,  decay  and  some  chemical  additive  effects  have  been  reported  in  the  1967 
Annual  Report/  A  few  additional  results  on  Lucite  are  reported  here 
together  with  new  results  on  thermoelectric  currents  in  PVC  and  in  a  few  other 
polymers.  Applied  field  and  EPR  results  are  presented  on  all  systems  studied. 

a.  Tharfreolostric  Currant  MmsuupiiqI 

Most  of  the  measurements  of  the  induced  thermoelectric  current  in  the 
organic  polymers  reported  here  were  made  on  disk-shaped  samples*  2.5  cm  in 
diameter  and  1.0  cm  thick.  The  samples  were  irradiated  to  a  level  of 
absorbed  dose  of  C.5  Mrad  and  subsequently  heated  at  a  constant  rate  of 
6°C/min.  Unless  stated  otherwise  these  values  should  be  assumed  for  all 
measurements. 

Thermoelectric  current  curves  for  Lucite  and  PVC  have  been  reported/ 
Figure  1-10  shows  the  thermoelectric  current  for  polystyrene.  The  curve  is 
composed  of  a  major  peak  of  1.1  pi  coamps  (pa)  at  1 1 5°C  and  a  plateau  of 
0.3  pa  between  60  and  80°C.  The  general  shape  of  this  curve  is  similar  to 
that  obtained  from  polystyrene  in  a  study  made  by  Hardtke  using  a  much 
larger  temperature  gradient. 


Figure  1-10.  Thermoelectric  current  plot  for  polystyrene. 

Doses  0.5  Mrad;  Heating  Rate!  6°C/minj 
Sample  Size:  2.0  cm  dia.  x  1.0  cm  thick. 
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Samples  of  polyethylene  and  of  Teflon  were  prepared,  irradiated  and 
heated.  No  measurable  thermoelectric  current  was  observed  in  repeated 
trials  of  either  of  these  materials  in  the  temperature  range  of  25  to  150°C. 

b.  Radiation  £oas  and  Thermal  Essay  Mnucimii  sn  EYC 

Samples  of  PVC,  1  cm  thick,  were  exposed  to  various  doses  of  Co60 
gamma-radiation  and  then  heated  in  a  temperature  gradient  to  produce  thermo¬ 
electric  currents.  Table  l-V  shows  the  results  of  measuring  the  currents 
generated  in  PVC  after  absorbing  doses  as  high  as  116  Mrad.  The  relationship 
between  normalized  total  charge  and  exposure  time  is  shown  graphically  in 
Figures  1-11  and  1-12  for  two  different  dose  ranges.  Values  of  normalized 
total  charge  for  exposure  times  to  1  hour  (0.5  Mrad)  are  shown  in  Figure  1-11 
while  data  for  exposure  times  to  22  hours  (11  Mrad)  are  illustrated  in 
Figure  1-12.  A  combination  of  the  curves  in  Figures  1-11  and  1-12  would 
consist  of  three  linear  segments  with  inflection  points  at  absorbed  doses  of 
approximately  0,12  Mrad  and  1.2  Mrad.  Saturation  of  the  trapped  charge 
carrier  population  does  not  occur  within  the  dose  range  studied. 

The  trapped  charges  in  PVC  are  quite  stable.  Samples  of  PVC  irradiated 
to  0.5  Mrad  and  stored  at  room  temperature  for  4  months  showed  no  observable 
decrease  in  the  normalized  total  charge  release  as  shown  in  Table  I -VI.  Data 
from  measurements  of  the  thermoelectric  current  in  PVC  after  the  release  of 
trapped  charge  carriers  at  various  temperatures  from  50  to  70°C  are  presented 
in  Table  t  —V II*  The  stability  of  the  charge  trap  appears  to  be  constant  to 
a  temperature  of  approximately  60°C  in  PVC.  This  is  illustrated  in  Figure  1-13 
where  the  normalized  total  charge  measured  in  PVC  samples  after  storage  for 
a  5  hr  period  at  various  temperatures  is  plotted  against  temperature. 


Table  l-V.  Effect  of  Absorbed  Dose  on  the  Thermoelectric 


Current  in  PVC* 


Exposuro 

Absorbed 

Peak 

Temperature 

Normalised  h 

Time 

Doso 

Current 

of  Peak 

Total  Charge0 

aea»— 

Mrad 

picoamps 

°C 

arbitrary  units 

5  min. 

.  04d 

5.8 

100 

14 

10  min. 

.03 

10.5 

102 

25 

15  min. 

.12 

13.5 

95 

32 

30  min. 

.25 

19 

97 

46 

60  min. 

,50 

31 

85 

67 

2  hr. 

1.0 

48 

90 

108 

3  hr. 

1.5 

70 

87 

134 

6  hr. 

3 

75 

83 

157 

12  hr. 

6 

100 

82 

199 

18  hr. 

9 

111 

80 

264 

22  hr. 

11 

134 

81 

270 

40  hr. 

20 

120 

84 

247 

78  hr. 

39 

152 

75 

349 

120  hr. 

60  , 

114 

75 

357° 

(two  peaks) 

64 

103 

186  hr. 

93 

120 

75 

367° 

(two  peaks) 

90 

93 

232  hr. 

116 

159 

78 

447° 

(two  peaks) 

129 

93 

^Sample  Size:  2.5  cm  dia.  X  1.0  cm  thick,  Heating  Rate;  6°C/*in 


b  o 

Curves  integrated  to  120  C 

cNormalized  Total  Charge  o t  both  peaks 

dUndorline  indicates  uncertainity 


EXPOSURE  TIME  (min) 


EXPOSURE  TIM 
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Table  l-VI.  Release  of  Trapped  Charge  Carriers 
in  PVC  at  25°C  a 


Decay 

Peak 

Tempo rature 

Normalized  ^ 

Time 

Current 

of  Peaks 

Total  Charge 

— — 

picoamps 

°C 

arbitrary  units 

10  hrs. 

28 

86 

56 

25  hrs. 

28 

87 

54 

4  days 

31 

SO 

61 

2  weeks 

27 

87 

60 

3  weeks 

28 

85 

58 

6  weeks 

26 

88 

57 

£  weeks 

25 

94 

55 

16  weeks 

27 

92 

61 

^Sample  Size: 

2,5  cn  dia. 

X  1.0  cm  thick,  Absorbed 

Dose;  0.5  Mrad 

h  C 

Curve  integrated  to  120  C 
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1-ihle  I -VII.  Release  of  Trapped  Charpe  Carriers 

a 

in  PVC  at  Various  1 emperatures 


Do  cay 

Time 

Peak 

Current 

Temperature 
of  Peak 

Normalized 
Total  Charge0 

hours 

pi  coamps 

°C 

arbitrary 

Temperature 

~50nC 

0 

31 

85 

67 

1 

30 

87 

66 

3 

37 

02 

64 

5 

37 

95 

62 

10 

35 

90 

54 

Temperature  60°C 


1 

3 

5 

7.5 

10 

43 

42 

42 

47 

43 

90 

92 

95 

93 

97 

67 

61 

62 

62 

63 

Tomporature 

6  5QC 

1 

50 

92 

70 

3 

43 

92 

58 

5 

37 

93 

48 

7.5 

33 

93 

47 

10 

28 

100 

35 

Temperature 

70“C 

1 

28 

96 

37 

3.6 

9.5 

99 

13 

5 

5.5 

95 

7 

*Sample  Size:  2.f  cm  dia.  X  1.0  cm  thick,  Absorbed  Dose:  0.5  Mrad, 
Keating  Rate:  6  C/min 


b  o 

Curves  integrated  to  120  C 


TOTAL  CHARGE  (arb.  units) 
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TEMPERATURE  (*C) 

Figure  1-13.  Five  hour  release  of  trapped  charge  carr:ers  at 
various  temperatures  in  PyC.  Dose:  0,5  Mradf 
Heating  Rate:  6°C/min$  Sample  Size:  2.5  cm 
dia.  x  1.0  cm  thick. 


Virtually  no  release  is  observed  after  5  hours  at  temperatures  below  60°C. 
However,  above  60°C  the  release  rate  increases  very  rapidly.  Complete 
curves  of  trapped  charge  carrier  release  for  temperatures  from  50  to  75°C 
are  shown  in  Figure  I —14- 

The  release  of  trapped  charge  carriers  in  PVC  follows  a  first-order 
rate  equation.  If  the  slopes  of  the  curves  of  trapped  charge  carrier 
release  at  60°C  and  above  are  determined  and  used  to  construct  an  Arrhenius 
plot,  Figure  1-15,  an  activation  energy  can  be  determined  from  the  slope  of 
the  straight  line.  A  value  of  67  kcal/mole  U.9  ev)  is  obtained. 

c.  Effect  fll  Blocking  Electrodes  fin  Thermoelectric  CudflOl 

Glass  blocking  electrodes  were  used  to  distinguish  the  flow  of  charge 
between  the  sample  and  the  external  circuit  from  the  internal  movement  of 
charge  carriers  within  the  sample,  lucite  disks,  1  cm  thick,  were 
irradiated  and  the  thermoelectric  currents  generated  in  the  disks  were 
measured  with  and  without  the  presence  of  blocking  electrodes.  Curves  of 
the  same  shape  and,  within  experimental  error,  of  the  same  magnitude  were 
observed  in  samples  measured  with  the  blocking  electrodes  as  was  observed 
in  samples  measured  without  the  blocking  electrodes.  This  is  evidence  that 
the  thermoelectric  current  is  due  entirely  to  a  recombination  of  charge 
carriers  within  the  volume  of  the  sample. 

d.  Thermoelectric  Currents  ifi  with  Addi.ti.YflS 

Additives  which  could  act  as  electron  or  hole  traps  were  doped  into 
synthesized  PMMA  to  study  their  effects  on  the  thermoelectric  currents. 
Additive  effects  on  photopolymer i zed  PMMA  have  been  reported  previously.^ 
The  effect  of  several  electron  scavengers  on  the  thermoelectric  current 


?0  2 
J08 


aoyvHo  ivioi 
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measured  in  thermally-polymerized  PMMA  are  shown  in  Table  I —V III.  The 
addition  of  biphenyl,  C02,  02  or  azobisisobutyronitrile  (AIBN)  produces  a 
significant  increase  in  thermoelectric  current  measured  in  the  synthesized 
polymer.  A I  BN  is  a  compound  which  also  functions  as  an  initiator  of 
polymerization  in  methyl  methacrylate.  Benzoyl  peroxide  added  before 
polymerization  likewise  produces  an  increase  in  the  thermoelectric  current. 
Since  benzoyl  peroxide  is  a  chemical  initiator  of  polymerization  commonly 
used  in  the  commercial  production  of  Lucite  it  must  be  at  least  one  of  the 
components  responsible  for  the  trapping  of  electrons  and  the  subsequent 
generation  of  thermoelectric  current  in  this  commerical  material. 

The  relative  effects  of  the  concentration  of  CCl4  on  the  generation  of 
charge  carriers  in  a  temperature  gradient  are  shown  in  Table  I  —  I X.  All 
of  these  samples  were  thermally  polymerized  and  then  machined  into  disks 
2.5  cm  in  diameter  and  0.4  cm  thick.  The  disks  were  irradiated  to  a  dose 
of  0.5  Mrad  prior  to  heating.  The  normalized  total  charge  generated  in  the 
samples  is  plotted  against  CCl^  concentration  in  Figure  1-16.  The  initial 
large  positive  slope  of  the  curve  indicates  that  a  very  small  concentration 
of  CCl4  has  a  large  effect  on  the  concentration  of  trapped  charge  carriers. 
The  normal izeo  total  charge  shows  a  linear  dependence  on  CCl4  concentration 
in  the  range  of  2  mole  %  and  8  mole  %  but  decreases  significantly  at  a 
concentration  between  8  and  12  mole  %.  The  decrease  at  12  %  may  be  due 
partially  to  extraneous  effects  of  the  polymerization  process  since  the 
scavenger  concentration  is  so  high. 

Disks  of  PMMA,  1.0  cm  in  diameter  and  0.4  cm  thick,  were  prepar&H  by 
thermally  polymerizing  MMA  with  various  additives  in  conce  trations  of 
1  mole  %  and  4  mole  %.  The  samples  were  irradiated  to  a  dose  nf  0.5  Mrad 
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Table  l-VII  I.  Effect  of  Electron  Scavengers  on  Thermoelectric 

a 

Current  in  Thermally  Polymerized  PMMA 


Additive 

Sample 

Thickness 

Peal: 

Current 

Temperature 
at  Peak 

Nor;  .all  zed  ^ 
Total  Charge 

(concentration) 

CM 

picoamps 

°C 

arbitrary  units 

Biphenyl 
(1  ir.ole  %) 

1.0 

5.6 

60 

161 

0 

1.0 

5.6 

55 

163 

(saturated) 

1.0 

5.9 

57 

174 

C02 

(saturated) 

1.0 

9.3 

63 

292 

AIBN 

0.4 

5.6 

63 

192 

(saturated) 

(two  peaks) 

2.6 

62 

cci 

(2  mole  %) 

1.0 

3.0 

55 

92.1 

0 

Sample  Size: 

2.5  cm  dia. , 

Dose:  o.5  Mrad,  Heating  Rate: 

6°C/min 

^Curves  integrated  to  120°C 


Table  IX.  Effect  of  CCl4  Concentration  on  Thermoelectric 


Currents 

£ 

in  Thermally  Polymerized  PMMA 

Peak 

Tempo raturo 

Normalized 

Average 

CC1.  cone. 

Current 

At  Peak 

Total  Charge** 

Normalized 

mole  % 

plcoamps 

°C 

arbitrary  units 

Total  Charge 

0 

1.5 

66 

38.3 

0 

1.1 

57 

24.9 

34,1  +  6.1 

0 

1.5 

60 

39.0 

a 

3.0 

55 

92.1 

92.1 

4 

3.-1 

59 

96.5 

91.9  +  4.6 

4 

3.1 

60 

87.3 

8 

3.8 

65 

105.8 

109.9  +4.1 

8 

3.9 

60 

114.1 

12 

1.6 

5S 

56.5 

55.4  +  1.1 

12 

2.1 

61 

54.3 

TOTAL  CHARGE  (arb.  units) 


ed  total  charge  vs.  CCl* 

.  Dose:  0.5  Mrad;  Heating 
e  Size:  2.5  cm  dia.  x  1.0 


and  then  heated  to  generate  thermoelectric  currents.  The  results  of  measuring 
these  currents  are  shown  in  Table  I-X-.  Increasing  the  concentration  of 
CCl^..  biphenyl,  and  2-methylpentene-1  by  a  factor  of  4  increases  the  average 
normalized  total  charge  measured  in  each  respective  polymer  by  a  factor  of 
approximately  two.  The  values  reported  in  Table  l-X  for  measurements  on 
samples  of  PMMA  containing  CCl^  are  smaller  than  values  reported  for  the  same 
materials  in  Table  I-  IXv.  This  is  because  of  the  correspondingly  smaller 
diameter  of  the  respective  samples.  -The  smaller  diameter  was  used  because 
cylindrical  sections  without  voids  could  be  produced  more  consistently 
during  the  polymerization  process. 

e.  Thermal ly-flenerated  Currents  in  an  Applied  Eififitcic  Field 

An  investigation  was  made  to  determine  how  an  electric  field  (E)  would 
influence  the  charge  carriers  generated  as  a  result  of  heating  a  previously 
irradiated  organic  polymer.  The  heating  of  all  samples  reported  in  this 
section  was  carried  out  in  a  zero  temperature  gradient  with  the  exception 
of  results  on  the  combined  effects  of  a  temperature  gradient  and  an  electric 
field. 

(1)  Electric  Field  Dependence 

The  results  of  apolying  a  potential  across  a  1.0  cm  thick 
irradiated  sample  of  Lucite,  while  heating  the  sample  uniformly,  are  shown 
in  Figure  1-17.  The  data  for  each  curve  was  measured  in  a  different  sample 
after  the  sample  had  been  given  a  gamma  dose  of  0.5  Mrad.  Electric  fields 
from  0  to  500  volts/cm  are  shown.  Each  of  the  curves  contain  a  peak  in  the 
temperature  range  of  100  to  115°C.  This  corresponds  approximately  with 

■3 

the  glass  transition  temperature  for  Lucite. ^  The  curves  do  not  return  to 
zero  current  but  rather  increase  sharply  above  120°C. 


Table  l-X.  Effect  of  Various  Additives  at  Different  Concentration* 
on  Thermoelectric  Currents  in  Thermally  Polymerized  PMMA 


l'cak 

f,'  .  ,  .  .  •  ,  ,  . 

i  O 

’.'or...  1. /<:>'. 

L> 

Total  Charge 

Additive 

cone 

Cll  *  i'C  v 

Av 

°c 

mole  % 

picoai.ips 

arbitrary  uni 

None 

0.2 

57 

6.7 

0.3 

52 

8.6 

CC1. 

4 

1 

0.7 

61 

19.9 

4 

1.7 

60 

51.2 

4 

1.5 

62 

40.9 

biphenyl 

1 

0.7 

60 

IS. 7 

1 

0.9 

58 

25.6 

4 

1.6 

60 

47.2 

4 

1.5 

62 

2-MP-l 

1 

0.5 

61 

13.7 

1 

0.6 

62 

15.3 

4 

.9 

65 

25.8 

4 

1.2 

68 

30.5 

CC1 -1  mole  % 

' _ t 


4  T 
and 

2-UP-l  1 

1.4 

64 

38.0 

1 

1.3 

5S 

35.6 

1.2 

60 

28.2 

4 

1.1 

62 

36.6 

*S ample  Sizo:  1.0  cm  dia,  X  1.0  cm  thick,  Dose;  0.5  Urad,  Heating  Rate: 
6&C/min 


b  o 

All  curves  integrated  to  120  C 


The  measurement  of  background  current  in  an  unirradiated  Lucite 
sample  during  heating  in  an  electric  field  is  shown  in  Figure  1-18  for  electric 
fields  of  90  and  500  volts/cm.  This  background  current  shows  a  significant 
increase  at  temperatures  above  90°C.  Background  current  measurements  made 
for  values  of  E  between  90  and  500  V  are  intermediate  to  the  curves  for 
these  two  values. 

In  order  to  separate  the  radiation-produced  affect  from  the 

background  current,  it  is  necessary  to  subtract  the  background  current  from 

the  curves  shown  in  Figure  1-17.  This  is  illustrated  for  the  90  and 

500  V/cm  curves  in  Figure  1-18.  The  result  of  this  correction  is  to  define 

the  peak  of  the  curves  better  and  to  give  them  a  more  uniform  shape*  The 
corrected  curves  are  shown  in  Figure  i-19.  The  temperature  at  the  peak  after 
the  correction  for  background  current  has  been  made  is  approximately 

100°C. 

The  current  peak  height  is  presented  as  a  function  of  the  applied 
electric  field  (E)  in  Figure  1-20.  There  is  a  linear  relation  between 
these  two  factors  for  values  of  E  to  1000  V/cm.  This  indicates  that  the 
current  is  ohmic. 

(2)  Radiation  Dose  Dependence 

Figure  1-21  presents  an  illustration  of  the  effect  of  accumulated 
exposure  to  the  gamma-flux  on  the  growth  of  charge  carrier  movement* 

Exposure  times  from  0  minutes  (background  current  only)  to  60  minutes  are 
shown.  This  corresponds  to  absorbed  doses  from  0  to  0.5  Mrad.  An  electric 
field  of  +300  V/cm  was  applied  for  each  measurement.  The  curves  are  not 
corrected  for  background  current.  The  integrated  current  versus  dose  data 
is  given  in  Figure  1-22  and  Table  I -X I  for  exposure  periods  up  to  4  hours 
(2  Mrad).  The  normalized  total  charge  reaches  approximately  95%  of 


E  =  500  V/CM 

CORRECTED 

CURVE 


/E*  90  V/CM 
CORRECTED 
j/CURVE 


E* 500 V/CM 
UNIRR  AD. 

SAMPLE  / 


E  >=90  V/CM 
UNIRR  AD. 
SAMPLE 


HEATING  TIME  (min) 


I 


CURRENT  (picoamps) 


HEATING  TIME  (min) 

Thermally  generated  current  profile*  corrected  for  background  current 
Doses  0,5  Mrad;  Heating  Rates  6°C/minj  Sample  Sizes  2.5  cm  dia.  x 
1.0  cm  thick;  Zero  Tmperature  Gradient. 


MAX.  CURRENT  (picoamps) 


Electric  field  vs.  maxi  stun  current  in 
Lucite.  Dose:  0.5  Mrad*  Heating  Rate: 
6®C/minj  Sample  Site:  2.5  cm  die.  x 
1.0  cm  thick*  Zero  Temperature  Gradient 


(VOLTS/ CM) 


HEATING  TIME  (min) 


Thermally  generated  current  profiles  in  a  300  V/ cm  field 
for  doses  below  saturation  in  Lucite.  Dose  Rate* 

0.5  Mrad/hrj  Heating  Rate*  6°C/min$  Sample  Size*  2.5  cm 
dia.  x  1.0  cm  thickj  Zero  Temperature  Gradient.  Curves 
labeled  for  exposure  time  in  minutes. 


OO  OOOOOO  QOOo 

srw  o  co  jo  j  w  Q  C0C05J-  oj 

(slmr^D)*"  39dVH0  “I  VI 01 
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saturation  at  a  dose  of  0.2f  Mrad.  This  is  to  be  compared  with  a  dose  of 
1.0  Mrad  for  the  saturation  of  normalized  total  charge  under  the  influence 


of  a  temperature  gradient. 

(3)  Polarization  Effect  in  Lucite 

During  successive  measurements  of  the  background  current  in  the 
same  electric  field  it  was  discovered  that  there  was  a  distinct  reproducible 
difference  in  this  current  if  the  electric  field  was  applied  across  the 
sample  while  the  disk  was  being  cooled  from  120°C  to  room  temperature. 

Figure  1-23  shows  the  results  of  five  successive  measurements  made  on  the 
seme  unirradiated  Lucite  sample  with  an  applied  field  of  500  V/ cm.  The 
only  difference  in  the  five  trials  is  whether  the  500  V  was  applied  during 
the  cooling  period  or  removed  at  the  temperature  of  120°C. 

Curve  A  in  Fig*  *  1-23  was  the  current  measured  during  the  first 
haatirg  of  the  sample.  The  electric  field  was  then  applied  continuously 
during  the  cooling  periods  after  trial  1  and  trial  2  and  curve  C  represents 
the  currents  measured  during  trial  2  and  trial  3.  At  the  upper  temperature 
of  trial  3 »  the  electric  field  was  removed  and  then  reapplied  to  the  sample 
after  it  had  reached  room  temperature.  Curve  B  represents  the  current 
measured  during  trial  4»  The  field  warn  again  removed  at  the  end  of  the 
heating  cycle  and  reapplied  at  room  temperature.  Trial  5  also  produced 
curve  B. 

The  results  of  the  five  current  measurements  represent  evidence 
that  the  background  current  measured  in  the  Lucite  samples  in  an  electric 
field  is  lower  if  the  electric  field  is  applied  during  the  cooling  period 
after  the  sample  is  heated  to  120°C.  This  suggests  that  the  background 


Figura  l-23» 


HEATING  TIME  (min) 

Effect  of  aiaetric  fiald  applied  during  cooling  on  background 
currant  in  unirradiatad  Lucita.  Firat  Haating*  Curval- 
f  appliad  during  cooling,  Second  Haating:  Curva  C-£  appl  ad 
duSSg  cooling.  Third  Haating:  Curva  C~€ 

cooling,  Fourth  Haating:  Ourva  B-€  flfil  «/«: 

cooling,  Fifth  Haating:  Curva  B.  Applied  Raid  (E):  jOO  V/«f 
Haating  Rata:  6°C/minj  Saeple  Sira:  2.5  cm  dia.  x  1.0  cm  thick. 
Zaro  Tamparatura  Grad i ant. 


current  may  be  due  to  some  type  of  a  polarization  effect  produced  by  the 

electric  field  during  the  heating  cycle.  Thie  effect  is  similar  to  that  used 
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in  the  production  of  electrets. 

(4)  Combined  Effect  of  Temperature  Gradient  and  Electric  Field 
The  effect  of  applying  both  the  temperature  gradient  and  an 
electric  field  to  a  sample  during  the  heating  cycle  was  investigated.  One 
cm  thick  Lucite  samples  were  irradiated  to  a  dose  of  0.5  Mrad  before  the 
measurements  were  made.  Figure  1-24  shows  the  results  of  combining  a 
♦200  V/cm  field  with  the  temperature  gradient.  Curve  A  represents  the  current 
in  a  Lucite  sample  while  under  the  influence  of  the  temperature  gradient 
only.  Curve  B  represents  The  current  in  a  similar  sample  under  the 
influence  of  a  +200  V/cm  field.  Curve  C  is  the  mathematical  sum  of  curves 
A  and  B  while  curve  0  is  the  current  measured  in  a  Lucite  sample  under  the 
influence  of  both  the  temperature  gradient  and  a  +200  V/cm  field.  The 
difference  between  curve  D  and  curve  C  is  probably  due  to  a  small  change  in 
the  temperature  gradient  produced  by  the  required  change  in  wire  connections 
made  to  the  sample  holder. 

A  measurement  was  also  made  with  the  temperature  gradient  and  a 
reverse  electric  field.  It  was  first  verified  that  a  reverse  electric 
field  would  produce  equal  and  opposite  currents  in  the  sample.  Nearly 
identical  curves  of  current  versus  heating  time  were  obtained  for  +  and  -300 
V/cm  fields  with  the  recorder  input  polarity  reversed  for  the  negative 
electric  field. 

Figure  1-25  shows  the  results  of  combining  the  temperature 
gradient  with  a  -300  V/cm  field.  Again  the  broken  curve  is  the  mathematical 
sum  of  the  two  individual  current  p’^ts  and  the  solid  curve  is  the  measured 
current  plot. 


CURREN  T  (picoampsj 


50 


3.0 


2.0 


30 


HEATING  TIME  fain) 


Figure  1-24. 


Combined  effects  of  temperature  gradient  and  electric  field  on 
current  generated  in  irradiated  Lucite.  Curve  Ai  Temperature 

gradient  only,  Curve  B:  E  field  only.  Curve  C:  Sum  of  A  and  B, 
urve  D:  Measured  curve,  E  field:  +200  V/cmj  Heating  Rate: 
6°C/minj  Dose:  0.5  Mradj  Sample  Size:  2.5  cm  dia.  x  1.0  cm  thick. 


CURRENT  (pic  o  amps) 


HEATING 


Combined  effects  of  temperaturf  gradient  and # reverse  electric 
field  on  current  generated  in  irradiated  Lucite.  —  Curve: 
Sum  of  individual  effect*,  -  Curve:  Measured  curve,  E  field 
-300  V/cm;  Heating  Rate:  6°Q/min}  Dose:  0,5  Mradj  Sample 
Sue:  2.5  cm  die*  x  1.0  cm  thick. 


70 


The  obvious  conclusion  from  Figures  1-24  end  1-25  is  that  the  two 
effects  are  additive.  Therefore  current  conduction  in  a  temperature  gradient 
and  in  an  electric  field  must  involve  different  conduction  mechanisms. 

(5)  Currents  from  Other  Polymers 

Several  polymers  were  observed  to  conduct  currents  in  an  electric 
field  but  not  under  the  influence  of  the  temperature  gradient  in  the  same 
temperature  range.  Figure  1-26  shows  current  measurements  made  while 
heating  pure  synthesized  PMMA  and  linear  polyethylene  (PE)  in  an  electric 
field.  Both  samples  were  1.0  cm  thick  and  were  given  a  dose  of  0.5  Mrad 
prior  to  measurement.  A  500  V/ cm  field  produces  about  the  same  current  peak 
height  in  pure  PMMA  as  a  300  V/cm  field  produces  in  PE.  The  current  peak 
from  pure  PMMA  occurs  at  66°C  which  is  approximately  the  same  temperature 
at  which  the  thermoelectric  current  peak  occurs  in  Lucite.  The  current 
peak  height  in  pure  PMMA  is  almost  twice  the  current  peak  height  measured  in 
Lucite  for  the  same  applied  electric  field.  The  peak  in  Lucite,  however, 
occurs  at  a  much  higher  temperature  than  the  peak  in  pure  PMMA.  The 
temperature  at  which  the  current  peak  occurs  in  PE  is  also  much  higher  than 
the  temperature  at  the  current  peak  in  pure  PMMA.  The  current  peak  in  PE 
occurs  at  a  temperature  of  90°C  which  is  20  to  40°C  below  the  melting 
temperature^  (Tm)  of  110  to  130°C  for  PE. 

Figure  I  —27  shows  the  current  conduction  in  Teflon  while  being 
heated  in  a  300  V/ cm  field.  The  sample  was  1.0  cm  thick  and  had  been  given 
a  dose  of  0.5  Mrad  prior  to  measurement.  A  peak  current  of  86  pi  coamps  at  a 
temperature  of  137°C  was  observed.  The  temperature  at  the  peak  is  well 
below  the  Tm  of  327°C  for  Teflon.  A  second  heating  of  the  same  sample  re¬ 
sults  in  the  background  current  shown  in  the  lower  curve  of  Figure  1-27. 
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Figure  1-26.  Thermally  generated  current  profiles  of  pure  PMMA  and 
PE  in  an  electric  field.  Doses  0.5  Mrad;  Heating  Rate 
6°C/min;  Sample  Sizes  2.5  cm  die.  x  1.0  cm  thick; 

Zero  Temperature  Gradient. 
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Nona  of  these  mater  ials—pure  PMMA,  linear  PE  or  Teflon— produced 
measurable  thermoelectric  currents  in  the  temperature  range  of  25  to  120°C. 

Irradiated  and  unirradiated  samples  of  PVC  were  also  heated  in  a 
zero  temperature  gradient  with  an  electric  field  applied.  The  current 
increased  continually  from  60°C  with  temperature  for  both  the  irradiated  and 
unirradiated  samples.  The  current  reached  the  maximum  range  of  the 
electrometer  (3  x  10”9  amps)  before  the  heating  cycle  was  completed.  Any 
radiation  effects  on  the  thermally-induced  currents  were  therefore  masked 
by  the  background  current. 


f.  MaftiU£flDfloi&  s£  Radical  fiacay  mXh  Temperature 

Samples  of  Lucite  and  PVC  were  irradiated  to  a  dose  of  0.5  Mrad  and 
their  free  radical  concentration  was  measured  by  EPR.  The  variation  of 
free  radical  concentration  was  observed  at  various  temperatures  while  the 
sample  was  heated  at  a  linear  rate  of  approximately  6QC/min.  Figure  1-28 
shows  the  results  of  such  a  measurement  on  the  two  different  sanples  of 
Lucite  (PMMA)  and  a  sample  of  PVC.  In  each  case  the  concentration  of  free 
radicals  begins  to  decrease  sharply  at  a  specific  temperature.  For  the 
two  Lucite  samples  this  temperature  was  approximately  65  and  720C, 
respectively*  and  for  the  PVC  sample  it  was  approximately  60°C.  The  initial 
increase  in  the  signal  intensity  from  the  Lucite  samples  during  the  heating 
process  was  unexpected  but  was  not  investigated  further. 


4.  DISCUSSION 

••  Modal  af  Ihsriroolgctrlfi  Current  Generation 
(l)  Charge  Mobility  and  Charge  Trapping 

The  results  from  the  synthesized  polymers  are  quite  revealing. 
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Figure  1-28.  Variation  of  EPR  signal  intensity  with 

temperature  in  gamma-irradiated  PMMA  and  PVC. 
PMMA  (Lucite),  Dose:  0.5  Mradj  Heating  Rate* 


which  contribute  to  the  thermoelectric  current,  appear  to  be  trapped  in  it 
during  irradiation  at  32°C.  Thus  either  or  both  positive  and  negative 
charges  are  mobile  during  irradiation  and  undergo  recombi nation.  If  a  good 
electron  trapping  agent,  such  as  CCl^,  is  present  the  thermoelectric  current 
increases  significantly.  This  suggests  that  at  room  temperature  in  pure 
PMMA  without  additives,  radiation-produced  electrons  are  mobile  and 
recombine  with  trapped  holes  during  irradiation.  If  the  positive  charge 
carriers  were  also  mobile  in  the  presence  of  an  electron  trap,  recombination 
could  still  occur,  and  no  increase  in  the  thermoelectric  current  would  be 
observed*  This  is  contrary  to  our  findings.  This  conclusion  is  also 
supported  by  the  experiments  with  2-methyl  pent  ene-1  (2-MP-1)  in  PMMA.  In 
solid  organic  matrices  2-4P-1  is  a  good  positive  center  trapping  agent  and 
a  relatively  poor  electron  trapping  agent.^  PMMA  doped  with  2-MP-1 
exhibits  a  relatively  small  current  increase  compared  to  pure  PMMA.  This 
result  would  be  expected  if  the  positive  centers  are  assumed  to  be  immobile 
during  irradiation. 

If  CCl^  or  benzoyl  peroxide  is  present  in  PMMA,  the  observed 
thermoelectric  current  increases  by  a  factor  of  approximately  10  over  that 
in  pure  PMMA.  CCl^  is  well  known  to  react  with  electrons  by  dissociative 
electron  attachment  in  solid  organic  matrices.  It  seems  likely  that  benzoyl 
peroxide  can  act  in  ix  similar  manner  because  of  its  weak  oxygen-oxygen 
bonds  and  the  electron  affinity  of  oxygen.  Alternatively,  it  is  probable 
that  peroxide  and  oxide  radicals  from  benzoyl  peroxide  are  produced 
during  polymerization  and  can  act  as  electron  traps.  It  should  be  pointed 
out  that  electron  traps  produced  by  dissociative  attachment  or  by  oxygen- 
containing  radicals  are  quite  thermally  stable,  and  may  well  be  stable  above 
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room  temperature  in  a  polymer  matrix.  The  results  reported  here  bear  this 
out.  Both  benzene  and  chlorobenzene  act  as  poor  electron  and  positive 
center  traps  in  solid  organic  matrices.  That  they  cause  no  increase  in  the 
thermoelectric  current  in  PMMA  is  also  consistent  with  the  results  cited 
above. 

Commercial  PMMA  (Lucite)  is  thermally  polymerized  by  using  benzoyl 
peroxide  as  a  chemical  initiator;  thus,  it  should  be  similar  to  PMMA  doped 
with  benzoyl  peroxide.  This  has  been  verified.  Lucite  is,  therefore, 
convenient  to  work  with  since  it  has  a  built-in  electron  trap.  It  should  be 
remembered,  however,  that  the  currents  observed  in  Lucite  are  characteristic 
of  added  electron  traps  and  not  of  the  pure  PMMA  polymer  itself. 

Polyvinyl  chloride  shows  a  relatively  large  current.  It  is  a 
polyalkyl  chloride  and  can  itself  react  with  electrons  by  dissociative 
electron  attacnment.  Alkyl  chlorides  readily  show  this  reaction  in  the 

n 

solid  phase.  The  large  thermoelectric  currents  observed  in  PVC,  compared 
with  PMMA  containing  2  mole  %  CCl^,  are  due  to  the  fact  that  the  number 
of  trapping  sites  is  larger.  Also,  the  C-Cl  bond  in  PVC  is  distributed 
uniformly  while  CCl^  in  PMMA  may  not  be.  Commercial  PVC  may  contain  addi¬ 
tives  which  also  act  as  electron  traps  but  their  effect  should  not  change  the 
behavior  of  PVC  itself. 

In  contrast  to  PVC,  Teflon  shows  a  negligible  thermoelectric 
current  in  the  temperature  range  studied.  There  is  no  trapping  in  this 
polymer  because  the  C-f  bond  is  strong  enough  to  preclude  dissociative 
electron  attachment. 

The  absence  of  a  thermoelectric  current  in  polyethylene  (PE)  above 
room  temperature  is  expected  because  H“  formation  is  improbable.  However, 
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Frankevich  and  Tal'roze  report  "bursts  of  emf"  in  electron  irradiated  PE 

plates  at  temperatures  both  below  and  above  room  temperature  while  the 

samples  were  heated  in  a  manner  which  resulted  in  a  temperature  gradient. 

Further  wr~k  on  PE  below  room  temperature  would  be  of  interest. 

(2)  Sign  of  the  Mobile  Charge  Carrier  * 

The  observation  r?  thermally  generated  currents  indicates  that 

one  or  both  charge  Carrie's  are  mobilized  by  thermal  stimulation.  Since 

the  currents  are  the  result  of  charge  carrier  movement  in  a  temperature 

gradient  they  are  thermoelectric,  and  the  most  mobile  charge  carrier 

9 

corresponds  to  the  polarity  of  the  colder  sample  electrode.  This  was  posi¬ 
tive  for  all  measurements.  Additional  evidence  that  the  positive  centers 
are  the  mobile  charge  carrier  is  provided  by  the  information  obtained  from 
the  use  of  additives  in  PMMA.  If  the  electrons  were  thermally  detrapped, 
the  temperature  at  which  the  peak  current  occurs  should  depend  on  the 
particular  electron  trapping  agent.  However,  the  temperature  at  which  the 
peak  current  occurs  is  about  the  same  for  all  electron  trapping  agents 
observed.  Therefore  it  must  be  the  trapped  positive  centers  and  not  the 


trapped  electrons  which  are  mobilized  by  thermal  excitation. 

Positive  ion  mobility  has  been  proposed  by  several  authors 
to  explain  their  measurements  of  conductivity  in  an  applied  electric 
field  during  irradiation  in  PE,  polystyrene,  Teflon  and  PVC.  Other 


10,11,12 


authors 
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propose  that  the  rad i at  ion- induced  conductivity  in  the  same 


systems  is  due  to  the  mobility  of  electrons.  These  assumptions  are  based 
on  an  attempt  to  explain  the  dependence  of  the  induced  current  on  the  dose 
rate  of  the  incident  radiation.  But,  none  of  these  authors  have  presented 
positive  experimental  evidence  for  the  identity  or  sign  of  the  mobile 
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charge  carrier  during  irradiation.  The  use  of  a  temperature  gradient  to 
generate  thermoelectric  currents  has  allowed  us  to  obtain  the  first 
unambiguous  sign  determination  for  charge  carriers  in  irradiated  pol^neric 
systems. 

b.  Charge  Trapping  Efficiency 

(1)  Estimation  of  G  Value  in  Lucite 

The  G  value  of  trapped  charge  carriers  in  Lucite  can  be  estimated 

from  the  saturation  of  normalized  total  charge  versus  dose  reported 

previously.1  The  normal! red  total  charge  generated  in  a  temperature 

gradient  saturates  at  a  dcse  ot  approximately  1.0  Mrad.  This  value  is  in 

agreement  with  the  dose  required  for  the  saturation  of  color  centers  in 

17 

PMMA  as  measured  by  Gardner  and  Epstein.  At  this  dose  it  can  be  assumed 
that  all  trapping  sites  in  the  material  have  been  filled.  It  has  been 
established  that  benzoyl  peroxide  provides  trapping  sites  in  Lucite. 

Benzoyl  peroxide  is  present  in  Lucite  in  a  concentration  of  eoproximately 
0.12  mole  %.  From  the  equation, 

*|qq  =  number  of  trapped  charge  carriers, 

where  D  is  the  dose  delivered  to  the  sample  in  ev,  we  can  calculate  G  if  it 
is  assumed  that  each  benzoyl  peroxide  molecule  acts  as  a  trapping  site  for 
one  charge  carrier.  The  dose  delivered  to  a  disk  shaped  sample,  2.5  cm  in 
diameter  and  1.0  cm  thick,  at  the  saturation  dose  of  1  Mrad  is  3.6  x  1G20  ev, 
and  the  number  of  molecules  of  benzoyl  peroxide  in  the  same  sample  at  a 
concentration  of  0.12  mole  i»  is  4-1  x  1019  molecules.  The  G  value  for 
trapped  charges  is  calculated  to  be  G  =  11.3.  The  validity  of  this  result 
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will  be  discussed  below.  This  G  value  is  based  on  the  assumption  that  each 
benzoyl  peroxide  molecule  forms  one  trap.  If  the  molecule  forms  two  traps, 
the  G  value  would  be  double.  The  calculation  also  assumes  that  the  only 
trapping  agent  is  benzoyl  peroxide.  Other  unknown  additives  in  the 
commercial  material  may,  in  fact,  contribute  to  the  trapping.  This  would 
also  raise  the  G  value. 

(2)  Calculation  of  the  Net  Distance  of  Charge  Carrier  Movement 

The  experiments  with  blocking  electrodes  indicate  that  the  charge 

conduction  is  due  almost  entirely  to  the  internal  recombination  of  charge 

carriers.  if  a  charge,  q,  moves  a  distance  d  between  two  parallel  plates, 

it  will  contribute  to  the  externally  measured  charge  Q,  where  D  is  the 

18 

total  distance  between  the  two  plates.  If  charge  carriers  move  at  an 
angle  other  than  90°  to  the  plates,  only  the  component  of  motion  perpen¬ 
dicular  to  the  plates  would  be  registernd  »n  the  external  circuit.  The 
movement  of  charge  carriers  under  the  influence  of  a  temperature  gradient 
is  a  process  of  diffusion  and  therefore  involves  random  motion  with  a  net 
movement  in  the  direction  of  the  temperature  gradient.  It  is  this  net 
motion  of  charge  carriers  that  is  observed.  If  the  externally  measured 
charge  is  Q,  then  the  average  component,  d,  of  net  movement  of  charge 
carriers  perpendicular  to  the  measuring  electrodes  equals  ^D.  The  value 
of  Q  measured  in  a  1  cm  thick  Lucite  sample  under  the  influence  of  a 
temperature  gradient  after  an  absorbed  dose  of  1  Mrad  corresponds  to 
2.6  x  10"9  coulombs.  The  value  of  q  is  equal  to  the  number  of  trapped 

charges  <•$&>  times  1.6  x  10  19  coulombs  per  charge  and  D  equals  1.0  cm. 

o 

From  the  equation  above,  d  was  calculated  to  be  1.2  A.  This  distance  seems 
unrealistically  small,  since  it  would  represent  movement  of  less  than  the 
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internuclear  distance  between  neighboring  atoms  in  the  polymer.  One  might 
therefore  conclude  that  the  calculated  G  value  wa s  much  too  high,  since  q 
is  dependent  on  G,  and  d  equals  ^  .  The  postulate  that  one  charge  carrier 
is  trapped  per  molecule  is  probably  invalid. 

If,  for  example,  10#  to  1#  of  the  benzoyl  peroxide  molecules  act  as 

trapping  sites,  the  distance  of  charge  movement  would  be  in  the  more 

o 

reasonable  range  of  12  to  120  A.  One  might  therefore  conclude  that  the 
trapping  efficiency  of  benzoyl  peroxide  in  Lucite  is  less  than  unity. 

Similar  arguments  apply  to  other  electron  trapping  agents.  Since  the  mobile 
charge  carrier  is  positive,  the  relative  efficiency  of  the  various  electron 
trapping  agents  is  proportional  to  the  normalized  total  charge  measured  in 
polymer  samples  containing  equal  concentrations  of  each  additive  after 
exposure  to  the  same  amount  of  radiation.  The  assumption  has  also  been 
made  that  one  positive  charge  carrier  participates  in  the  thermoelectric 
current  for  every  trapped  electron.  If  only  a  fraction  of  the  trapped 
positive  charge  carriers  is  able  to  participate  in  the  thermoelectric  current, 
the  calculated  recombination  distance  would  have  a  higher  value. 

(3)  Comparison  of  Lucite  with  PVC 

Polyvinyl  chloride  (PVC)  gives  much  larger  currents  in  a  tempera¬ 
ture  gradient  than  does  doped  PMMA  primarily  because  of  the  larger  number 
of  trapping  sites  in  this  material.  The  normalized  total  charge  in  the 
sample  does  not  saturate  at  doses  as  large  as  100  Mrad.  The  normalized 
total  charge  measured  in  PVC  under  the  influence  of  a  temperature  gradient 
after  a  dose  of  0.5  Mrad  is  1.7  x  10”8  coulombs.  If  a  G  value  of  0.1  is 
assumed,  5.7  x  10”2  couloumbs  would  bo  trapped  in  the  sample,  and  each 
mobile  charge  carrier  would  move  an  average  of  30  %  before  recombination 
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occurred.  Again  this  indicates  that  a  value  of  G  much  less  than  1  must  be 
assumed  in  order  to  obtain  a  reasonable  value  of  d.  Since  G  for  total 
ionization  is  normally  around  3  we  see  that  only  a  small  fraction  of  the 
initially  formed  ions  are  trapped. 

c.  Model  Conduction  Mechanism 

3  .... 

Deutsch  and  co-workers  have  measured  the  dispersion  in  the  dielectric 

constant  and  the  mechanical  compliance  of  PMMA  at  various  temperatures.  By 
plotting  log  frequency  versus  temperature  curves  for  the  various  parameters 
they  obtain  the  loci  of  maximum  change  in  these  properties.  From  these 
studies  they  have  identified  two  transitions  in  PMMA  which  they  have  labeled 
alpha  and  beta.  Alpha  is  the  transition  at  the  highest  temperature,  120°C, 
and  corresponds  to  the  change  from  the  glassy  to  the  rubbery  state.  It  is 
caused  by  the  movement  of  large  groups  of  atoms  or  the  main  polymer  chain. 
Values  for  the  softening  point  or  glass  transition  temperature,  Tg,  in 
PMMA  are  quoted  by  different  authors  as  temperatures  ranging  from  90°C  to 
120°C.  This  range  of  temperatures  may  be  due  to  different  degrees  of 

polymerization  and  to  the  fact  that  different  physical  parameters  were 
measured  as  an  indication  of  the  transition  point. 

The  beta  transition  in  PMMA  was  observed  by  Deutsch,  et  al.,  at  50°C. 
This  was  determined  from  measurements  of  the  dielectric  loss  of  the  polymer 
at  various  temperatures.  The  change  in  dielectric  loss  is  ascribed  to  the 
onset  of  free  rotation  of  the  polar,  carbomethoxy  side  groups  attached  to 
the  main  polymer  chain.  From  plots  of  log  (frequency  at  loss  maximum) 
versus  l/T°K,  activation  energies  (E&)  were  calculated.  Values  for  the  beta 
transition  in  PMMA  were  determined  by  measurements  of  two  different  physical 


parameters.  From  dielectric  loss  measurements,  a  value  of  20  kcal/moie  was 
calculated,  and  from  audio  frequency  mechanical  loss  measurements,  was 
determined  to  be  18  kcal/mole.  This  is  in  excellent  agreement  with  our 
measured  activation  energy  of  18.6  kcal/mole  calculated  from  the  thermal 
decay  of  trapped  charge  carriers  in  PMMA.  Furthermore,  the  temperature  at 
which  the  peak  thermoelectric  current  occurs  in  Lucite  is  comparable  to  the 
temperature  of  the  beta  transition.  This  correlation  of  activation  energy 
and  peak  temperature  strongly  suggests  that  the  detrapping  of  charge  carriers 
in  PMMA  is  associated  with  the  rotation  of  the  carbomethoxy  side  groups. 

One  can  then  envision  an  attractive  conduction  mechanism  as  illustrated 
below  where  positive  centers  are  transferred  between  adjacent  side  groups  and 
thus  propagated  through  the  sample. 

CH3  ch3 

I  Y  I 


II  I  I  II 
o  ch3  h3c  0 


This  conduction  mechanism  should  show  stereospecificity  and  depend  on  whether 
adjacent  carbomethoxy  groups  are  on  the  same  side  (isotactic  polymer)  or 
opposite  sides  (syndiotactic  polymer)  of  the  main  chain.  The  polymer  we 
have  used  was  largely  atactic  (random  orientation  of  side  groups).  Future 
experiments  are  aimed  toward  testing  the  stereospecificity  of  the  thermo¬ 
electric  current  magnitudes  and  the  trapped  charge  stability  in  PMMA. 

There  is  also  a  correlation  between  the  temperature  of  the  alpha 
transition  in  PMMA,  between  90°C  and  120°C,  and  the  temperature  at  which 
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the  maximum  current  occurs  in  PMMA  in  an  applied  electric  field.  This 
suggests  that  the  population  of  rad i at  ion- induced  charge  carriers,  seen  with 
an  applied  electric  field,  are  not  mobile  until  a  temperature  is  reached  at 
which  the  main  polymer  chains  can  move. 

The  peak  thermoelectric  current  measured  in  PVC  also  occurs  approxi¬ 
mately  at  its  glass  transition  temperature.  In  this  case  there  is  no  lower 
temperature  transition  due  to  internal  motion  and  main  chain  motion  is 
apparently  necessary  for  conduction  to  occur.  We  have  found  no  activation 
energy  measurements  in  PVC  with  which  to  compare  our  values.  The  h.'gh 
temperature  stability  of  the  trapped  charge  carriers  in  PVC  together  with 
the  ability  to  build  up  trapped  charge  without  saturation  to  very  large 
doses  (to  at  least  100  Mrad),  suggests  that  the  PVC  might  have  possible 
application  where  it  is  desirable  to  store  relatively  large  amounts  of 
charge  for  long  periods  of  time. 

d.  Thermally  Generated  Charge  Conduction  in  an  Electric  Eiaid 

The  results  reported  in  Section  3e  show  conclusively  that  the  conduction 
mechanisms  and  mobile  charge  carrier  populations  are  different  for  currents 
thermally-generated  in  an  electric  field  compared  with  currents  produced  by 
a  temperature  gradient  in  irradiated  lucjte,  The  explanation  for  this 
difference  is  not  readily  apparent.  More  definitive  experiments  will  be 
necessary  before  an  adequate  description  of  the  mobile  charge  population 
responding  to  the  electric  field  can  be  given. 

The  effect  of  electron  scavengers  on  the  radiation-induced  electrical 
conductivity  should  be  studied  in  greater  detail.  How  does  the  concen¬ 
tration  of  these  additives  influence  the  current  and  the  temperature  at  which 
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the  maximum  current  occurs?  Why  is  a  current  observed  upon  heating 
irradiated  samples  of  some  materials,  such  as  PE,  Teflon  and  pure  PMMA, 
in  an  electric  field,  but  no  thermoelectric  current  is  observed  in  these 
materials  in  the  same  temperature  range?  These  and  other  questions  need  to 
be  answered. 

One  would  expect  an  electric  field  to  influence  the  movement  of  the 
same  thermally  liberated  charge  carriers  that  would  be  influenced  by  a 
temperature  gradient.  The  observation  that  the  currents  produced  by  these 
two  current  generating  mechanisms  are  independent,  and  can  be  added  or 
subtracted,  is  very  puzzling.  One  possible  explanation  for  the  origin  of 
the  current  produced  by  the  electric  field  is  that  the  radiation  polarizes 
soma  part  of  the  polymer  chain  which  remains  polarized  until  the  sample  is 
heated  above  the  glass  transition  temperature.  At  this  temperature 
polarization  relaxation  could  provide  the  observed  current.  Further  in¬ 
vestigation  of  these  phenomena  is  in  progress. 
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In  a  continuing  theoretical  study  concerned  with  the  characterization  of 
highly  energetic  species  having  potential  interest  to  the  U.  S.  Air  Force, 
the  tools  of  ab  initio  quantum  mechanics  have  been  used  both  to  characterize 
several  individual  molecules  and  to  devise  new  techniques  of  general 
applicability  to  molecular  problems.  Studies  of  this  nature,  although 
time-consuming  and  often  difficult,  provide  exceedingly  powerful  tools  for 
the  understanding  of  the  chemistry  of  small  molecules,  since  the  solutions 
of  the  Schrodinger  equation  apparently  provide  the  only  infallible  basis 
for  the  discussion  of  chemical  phenomena  associated  with  isolated  molecules 
or  small  groups  of  them. 

The  individual  species  studied  include  H3+,  HeH* ,  He2++,  H2  and  Li+. 
These  species  were  examined  from  the  point  of  view  of  establishing  a  lower 
bound  estimate  of  the  energy  of  the  ground  state.  This  is  an  especially 
important  study  for  molecules  such  as  H3+,  where  the  magnitude  of  the 
binding  energy  is  crucial  to  discussions  of  possible  stabilization  by  anions, 
and  is  not  known  from  experiment. 

The  effort  directed  toward  developing  new  and  innovative  techniques 
which  will  allow  quantum  chemists  to  perform  calculations  which  are 
accurate  enough  for  the  prediction  and  interpretation  of  phenomena  such 
as  chemical  bonding  has  been  carried  on  in  two  separate  directions. 

In  order  to  make  lower  bound  estimates  for  multi -electron,  polyatomic 
molecules,  it  was  necessary  to  devise  new  mathematical  techniques  for 


I 
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handling  the  problem.  It  wa9  found  that  Gaussian  basis  orbitals  provided  an 
excellent  answer  to  this  search,  for  all  of  the  integrals  that  arise  can 
be  evaluated  either  in  closed  form,  or  by  a  single  numerical  integration. 
Furthermore,  the  application  of  this  technique  provided  a  very  sensitive 
test  of  a  trial  wavefunction,  and  indicates  where  improvements  should  be 
made. 

The  development  and  implementation  of  tfor.  genunal  theory  of  chemical 
bonding  also  experienced  major  developments,  ihe  new  technique  for  the 
solution  of  the  coupled  integro-differentiM  equations  that  arise  in  the 
determination  of  the  optimum  geminals  for  a  system  has  been  developed 
considerably,  and  an  iterative  method  having  quadratic  convergence  to  the 
true  solution  has  been  developed.  In  addition,  the  connection  between 
geminal  theory  and  traditional  Hartree-fock  theory  has  been  explicitly 
derived,  so  that  a  direct  assessment  of  the  expected  improvement  in  accuracy 
can  be  made. 

The  implementation  of  the  calculations  using  geminal  theory  on  the 
GE-625  has  proceeded  slowly,  and  the  integral  programs  are  now  fully 
operational.  Since  small  inefficiencies  in  programs  of  this  nature  result 
in  large  increases  in  the  amount  of  computing  time  necessary  to  perform 
these  computations,  it  was  deemed  essential  to  proceed  slowly  in  their 
development,  to  be  certain  that  efficiency  had  not  been  sacrificed.  The 
consummation  of  these  efforts  in  terms  of  the  programs  to  determine  the 
energy  and  geminals  using  these  integral  values  is  now  underway. 


B.  LOWER  BOUND  AND  ENERGY  VARIANCE  CALCULATIONS 
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1 .  PURPOSE 

The  Variation  Theorem  and  methods  of  computation  such  as  configuration 
interaction  provide  powerful  techniques  for  obtaining  the  lowest  energy 
possible  for  a  wavef unction  with  adjustable  parameters.^  The  power  of 
these  techniques  lies  in  the  fact  that,  for  an  arbitrary  choice  of  the 
adjustable  parameters,  the  energy  associated  with  the  trial  wavefunction 
will  always  be  above  the  true  energy  of  the  system.  The  difficulty  in  such 
an  approach  is  that  the  method  in  which  the  trial  wavefunction  should  be 
found,  i.e.,  which  terms  and  how  many  terms  should  be  included,  is  difficult 
to  formulate  in  general.  Thus,  there  is  a  great  need  for  a  technique  that 
will  provide  more  information  concerning  the  accuracy  of  both  the  approxi¬ 
mate  energy  and  the  associated  wavefunction  at  any  stage  of  complexity. 

It  is  known  that  an  examination  of  the  energy  variance  associated  with 

a  trial  wavefunction  provides  a  means  by  which,  in  principle,  this  desired 

2 

information  can  be  obtained.  The  variance  is  defined  as 

U2  =  2>  -  >2  £  0,  (1) 

where 

<ti*>  =  j  <ww*)  *,  <2) 

and 

<#>  =  5  <p*^<p  (3) 

Here  <p  is  a  normalized  trial  wavefunction,  is  the  Hamiltonian  operator 
for  the  system  under  consideration,  and  the  integration  is  performed  over 
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toe  coordinates  of  all  of  the  electrons.  Because  of  the  nature  of  the 
operators  contained  in^f2,  the  energy  variance,  which  should  be  zero  for 
an  exact  solution,  provides  a  sensitive  test  for  any  proposed  trial  function. 
Also,  the  formulas  of  Temple,  Stevenson*  and  Weinstein  can  be  used  to 
estimate  a  lower  bound  to  the  true  energy  of  the  system  if  the  variance  is 
known.  Consequently,  an  evaluation  of  U2  provides  not  only  a  means  of 
estimating  the  overall  accuracy  of  the  trial  wavefunction,  but  also  a  means 
of  bracketing  the  true  energy  of  the  system. 

When  Slater-type  orbitals  are  used  as  basis  orbitals,  the  evaluation 
of  the  integrals  contained  in  <^f2>  is  exceedingly  difficult.^  If,  on  the 
other  hand,  the  trial  function  is  formed  using  a  basic  set  consisting  of 
Gaussian-type  orbitals,  the  necessary  integrals  over  the  square  of  the 
Hamiltonian  can  be  calculated  both  accurately  and  conveniently. 

In  this  study,  the  energy  variance  for  a  series  of  simple  systems  has 
been  examined  using  Gaussian  orbitals.  The  systems  that  were  investigated 
are  H2+,  H2,  HeH+,  He2++  and  H3\ 

2.  ENERGY  VARIANCE  CALCULATIONS 

The  wavefunctions  for  the  above  systems  that  were  employed  in  this 

7,8 

study  have  been  reported  recently  by  M.  Schwartz  and  L.  Schaad.  In  that 

study,  a  basis  set  consist ing  of  unnormalized,  Is-type  Gaussian  orbitals 

9 

was  employed, 

Ga(1)  =  exp(-ar2A1),  14) 

where 

r2Ai  =  (x,-Ax)2  +  (yj-Ay)2  +  lz,-Az)2 


9 


(5) 
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and  where  "a"  is  a  variable  orbital  exponent.  Also,  (Ax,  Ay,  A2)  locate  the 
position  of  the  Gaussian  orbital  relative  to  some  arbitrary  origin.  The 
coordinates  (Ax,  Ay,  Az)  do  not  necessarily  coincide  with  a  nucleus,  thus 
allowing  for  the  possibility  of  "floating"  orbitals. 

In  the  determination  of  the  wavefunction  to  be  used  here,  the  energy  was 
minimized  with  respect  to  all  the  parameters,  including  the  position  and  size 
of  the  Gaussians.  The  non-relativistic  Hamiltonian  for  a  syvtem  of  N 
electrons  and  M  nuclei,  within  the  Born-Oppenheimer  approximation,  can  be 
written  in  atomic  units^  as 


N  N  M  7  N  M  77 

i=i  i=i  A=i  A1  i<j  ,J  A<B 


(6) 


where  the  indices  i  and  j  refer  to  electrons,  and  A  and  B  refer  to  nuclei. 

By  direct  substitution  of  equation  (6)  into  equation  (2^,  it  can  be 

seen  that  the  new  integrals  that  arise  in  the  evaluation  of  <^f2>  that  have 
not  already  been  evaluated  in  <if>  are,  when  written  in  their  most  general 
form, 

J,  =  5t  V,2<pA(l)][V,2q>B(1)]  dv,  ,  (7) 

ja  =  str-^dlH^VU]  dv,  ,  (8) 

Ci 

J,  =  55[v,2  (pAd)  <pcU)][-~  <Pr(1)  <pd(2)]  dv,  dv2  ,  (9) 

»  n  r  j  2 

^  ~  JCtJ-  <Pa0)][~  <PB^)3  dvt  » 


(10) 
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o,  =  jj[t^  <pa(D  <pc(2)]P-  *B(1)  »du>]  dv>  dv>  •  l11> 

J4  =  J[^-<pA(1)][^-<P0O)]  dv,  ,  (12) 

J,  =  5$[^;<pA(i)  9c(2)][^Tb(1)  »D(2)]  dv,  dv,  ,  (13) 

J«  -  JiJ E-pi—  9A(1)  <Pq(2)  cpE(3) jc^r*—  <Fg(l)  dp (2)  <PF (3)]  dv,  dva  dv3  . 

(14) 

The  reason  for  the  intractability  of  these  integrals  when  Slater -type 
orbitals  are  employed  is  apparent,  due  to  the  presence  of  a  four-center 
one-electron  integral,  as  well  as  a  five-center  two-electron  and  a  six-center 
three-electron  integral.  A  discussion  of  the  methods  used  to  evaluate  these 
integrals  is  given  in  the  subsection  on  integral  evaluation. 

For  Ha  ,  the  nine-term  optimized  wavefunction  for  the  ground  state 
was  employed.  For  the  diatomic  molecules  H2,  He2  and  HeH  ,  the  wave- 
functions  as  given  in  Tables  IX,  XI  and  Xi I  of  Reference  (7)  were  employed. 
These  configurations  take  into  account  the  left-right  and  in-out  correlation 
of  the  electrons,  with  a  complete  optimization  of  all  the  parameters. 

Finally,  for  H3+  the  multi-conf igurational  wavefunct ions,  whose  complexity 
is  defined  by  the  parameters  (h.,  n2,  n3),  were  used.  These  calculations 
on  H3+  represent  the  best  variational  calculations  on  this  molecule  to  date. 
The  total  wavefunction  for  the  two-electron  molecules  can  be  written  as 

♦  (1,2)  =  <^(1)  C  <p(2)  H5) 

A#  <V  M 

where  £  is  a  column  vector  of  Is-type  Gaussian  orbitals,  and  the  coefficient 
matrix  fi  is  real  and  symmetric. 


The  results  of  the  computations  using  these  wavefunctions  are  giver  in 
Table  ll-l.  Included  in  this  table  are  the  expectation  values  over  and 
tf7,  as  well  as  the  variance  U2,  the  Weinstein  lower  bound  and  the 
expectation  values  of  the  three  components  of  2f2,  T2,  (TV+VT),  and  V2, 
respectively. 

3.  DISCUSSION 

The  value  of  the  energy  variance  is,  with  the  possible  exception  of  H2+, 
far  from  the  desired  zero  result.  This  is  not  particularity  surprising, 
since  the  parameters  contained  in  the  wavefunction  have  been  obtained  by  the 
minimum  energy  condition  of  the  variational  method.  This  does  not  neces¬ 
sarily  guarantee  that  the  expectation  value  of  the  other  operators  will  be 

predicted  accurately,  since  another  quantum  mechanical  operator  may  weight 

11 

a  different  region  in  space  than  does  the  Hamiltonian.  Thus,  it  is 

expected  that  the  operators  contained  in  «£f2 (=  T2  +  TV  +  VT  +  V2)  will 

provide  a  critical  test  of  any  wavefunction.  Operators  such  asV(4  depend 

not  only  on  the  accuracy  of  the  charge  density  close  to  the  nuclei,  but 

also  on  the  slope  of  the  wavefunction  in  that  region.  It  is  in  this  region, 

however,  where  Gaussian  orbitals  generally  have  difficulty  providing  a  good 

description  of  the  charge  distribution,  which  is  a  point  that  has  been 

12 

emphasized  recently  by  Schwartz.  In  the  case  of  the  hydrogen  atom, 
Schwartz  compared  the  local  energy  for  the  exact  solution  to  that  of  a 
Gaussian  approximation  to  this  exact  solution.  Although  tht>  energy  was 
within  10  hartrees  of  the  true  energy,  the  best  energy  variance  was  only 
0.013  hartrees2.  It  was  pointed  out  that,  for  the  hydrogen  atom  case  in 
particular, 
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12  13 

should  be  constant  for  an  exact  solution  *  and  that  deviations  from 
constancy  for  a  trial  function  can  be  attributed  to  a  variance  in  the  terms 
containing  the  momentum  operator.  Thus,  one  might  expect  for  the  more 
complicated  molecules  under  consideration  here  that  this  situation  will  be 
magnified.  In  addition,  any  inadequacy  of  the  wavefunction  to  describe  the 
electron  correlation  is  also  expected  to  be  a  contributing  factor  to  the 
magnitude  of  the  energy  variance. 

Unfortunately,  it  is  not  possible  to  compare  the  components  of  <tf2>} 
namely,  <T2>,  <TV+VT>  and  <Va>,  with  their  "true"  values  for  the  cases  under 
consideration  here.  However,  a  comparison  of  these  quantities  as  a  function 
of  the  complexity  of  the  wavefunction  can  give  some  indication  of  their 
convergence.  For  this  reason  they  are  also  listed  in  Table  ll-l.  As  can 
be  seen  for  those  cases  where  a  variety  of  trial  wavefunctions  were  used, 
the  rather  rapid  convergence  of  the  energy  is  not  accompanied  by  a  corres¬ 
ponding  convergence  of  the  energy  variance. 

The  investigations  that  were  made  on  the  one-electron  H2+  ion  provide 
an  energy  variance  of  0.023  hartrees2,  which  is  comparable  in  magnitude  to 
that  obtained  for  the  hydrogen  atom  (.013  hartrees2).  Although  the 
energy  in  both  these  systems  is  accurately  predicted,  the  magnitude  of  the 
quantity,  <ff2>  -  </f>2,  is  far  from  zero.  This  result  emphasizes 
dramatically  that  accentuates  the  features  of  the  wavefunction  not 

accounted  for  in  its  variational  determination.  Furthermo'e,  it  points  out 
how  the  supposition  that  a  good  energy  necessarily  leads  to  an  accurate 
overall  wavefunction  can  be  erroneous. 

The  calculations  on  the  two-electron  systems  Ha  and  H3+  provide  even 
more  striking  examples  of  the  inability  of  the  minimization  of  the  upper 
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bound  to  provide  a  good  energy  variance.  In  these  cases  a  decrease  in  the 
accuracy  of  the  total  energy,  due  primarily  to  the  difficulty  of  describing 
tho  electron  correlation,  is  reflected  in  a  marked  increase  in  the  energy 
variance  over  that  found  for  the  one-electron  systems.  Furthermore,  the 
sensitivity  of  U2  to  the  complexity  of  the  wavefunction  is  clearly 
demonstrated.  For  example,  a  consideration  of  the  two  best  wavefunctions 
for  the  hydrogen  molecule  shows  that  the  change  in  U2  is  100  times  greater 
than  the  change  in  <tf>.  Similarity,  for  H3+  a  change  in  <^f>  of  .001 
hartrees  is  accompanied  by  a  change  in  U2  of  .028  hartrees2.  Also,  the  large 
change  in  the  variance  cannot  be  traced  directly  to  any  one  of  the  components 
of  <2f2>,  for  they  are  all  far  from  constant,  although  <V2>  does  appear  to 
be  converging  slightly  faster  than  <TV+VT>  or  <T2>.  Thus,  the  extension  to 
the  case  of  two  electrons  with  the  concommitant  introduction  of  correlation 
effects  suggests  that  the  poor  lower  bound  is  not  only  due  to  the  lack  of 
accuracy  in  the  description  of  the  local  kinetic  energy  near  the  nucleus, 
but  also  that  <V2>  is  poorly  represented. 

A  similar  situation  exists  for  He2++  and  HeH+.  In  addition  to  the 
considerations  already  discussed,  the  larger  nuclear  charges  are  more 
demanding  of  the  trial  function  in  the  vicinity  of  these  nuclei.  The  inability 
of  the  Gaussian  description  to  accurately  reproduce  both  the  slope  and 
magnitude  of  the  wavefunction  in  this  region  is  responsible  for  the  electronic 
energy  of  these  systems  being  considerably  higher  than  their  true  values. 

Also,  it  is  these  features  of  the  Gaussians  that  are  heavily  weighted  by  the 
operator  J{2.  As  a  consequence,  the  components  of  <2f2>  are  greater  than 
those  of  Ha  by  an  order  of  magnitude,  and  their  lack  of  convergence  is  not 
unexpected. 
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To  estimate  the  accuracy  of  the  charge  distribution  around  a  nucleus, 
a  confutation  of  the  forces  provides  a  useful  test.  For  a  heteronuclear 
diatomic  molecule,  AB,  the  net  force  on  nucleus  A  in  the  z-direction,  equal 
to  -iE/azA,  is  given  by1^ 

Fa  =  FzAn  (nuclear)  +  FzA®  (electronic),  ^7) 


or 


according  to  the  coordinate  system  in  which  the  origin  is  chosen  at  atom  A, 
and  +z.  is  measured  in  the  direction  toward  atom  B.  Fg  is  similarity  defined. 

n 

Here  the  perpendicular  forces  on  the  nuclei  A  and  B  are  zero  by  symmetry. 

The  first  term  in  Equation  (18)  represents  the  nuclear  force  of  repulsion, 

and  the  second  term  represents  the  electronic  force  of  attraction.  For  H3 
the  situation  is  rather  more  complicated*  Here  the  total  force  can  be 
divided  into  two  components,  F|(  and  F^  respectively,  as  show  in  Figure  11-1. 

2 

F  *  -2-1  -  $  +  (1,2)  l  -ii?  +(1,2)  dx,  dx,  ,  .  (19) 

"  *rAB  rAi 

and 

2  y 

F  =  -f^r-J+0,2)  l  -2^y  +  (1,2)  dx,  dx,  .  (20) 

1  2rAB  i='  Ai 


For  a  stable  molecule  the  net  forces  on  the  nuclei  must  be  zero.  Thus 


fa  -  Fe  - 


F  =  F  =0. 
r„  r± 


(21) 


V 


i k£i 
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The  value  of  these  resultant  forces  for  the  approximate  wavefunction*  are 
listed  in  Table  1 1— 1 1 •  These  values  are  generally  very  good,  and  come  close 
to  meeting  this  requirement,  with  the  obvious  exception  of  HeH*,  and  to  a 
lesser  extent  He*4*.  In  HeH*  the  large  positive  value  of  F  ,  in  contrast 
to  the  nearly  zero  value  of  F^,  indicates  that  there  is  a  deficiency  of 
charge  to  the  bonding  side  of  the  hydrogen  nucleus.  This  suggests  that,  in 
the  determination  of  the  wavefunction,  the  larger  nuclear  charge  associated 
with  the  Helium  nucleus  places  a  more  stringent  condition  on  the  energy, 
which  the  wavefunction  tries  to  meet  at  the  expense  of  the  description  of 
the  charge  associated  with  the  proton*  Such  a  deficiency  of  charge  is 
undoubtedly  a  largo  contributing  factor  to  the  poor  value  of  the  energy 
variance.  In  He2**  the  higher  nuclear  charges  are  responsible  for  the  value 
of  F^  to  be  an  order  of  magnitude  higher  than  F^  in  Ha,  which  is  a  result 
that  is  consistent  vith  the  increase  in  magnitude  of  <#*>  with  increasing 
nuclear  charge. 

4.  INTEGRAL  EVALUATION 

The  integrals  Jf,  Ja  and  J,  can  be  readily  obtained  from  the  overlap, 
nuclear  attraction  and  electron  repulsion  integrals,  respectively,  by 
recalling  the  identity 

V,2  Ga(1)  =  ?A2  Ga(1)  (22) 

where  now^A2  is  independent  of  the  electronic  coordinates,  and  hence  the 
integration  variables.  Thus,  we  have 


.0*2*7 


Ji  -  7/  7b*  J  cA(l)  Gg(l)  <»vi  # 

■  *( ;) *  (f)*  £  «>■*>(  f )  ®*  ♦ 46  (f)*  «*  •  “*  l  • 

J.  -  V,*  I  6.(1)  7-  0,(1)  On  , 

■  Cl 

-  t  ?b-5^-F«(Xx)  +Q^)aKP  J?i(Xx) 

j>  -  V  JJ  0A(1)  «c(a)  ^  6,(1)  6,(2)  4.1  0»,  , 

t  8-b-SPQ  -  6  ]  Fx(Xi)  ♦  £ 


vbare 


p  ■  (a+b)  , 

?t  «  (aAt+b»i)/p  . 


and,  in  general, 


q  ■  (c+d)  , 

Qt  ■  (cCj+dD^J/q  , 


**»!»* 


»* -3;*  <v  V*  •  "■*S<V,x)(V,i) 


Alto,  K  =  exp  ~  •  ABa)  , 

P 


X,=  pCPa  , 


L  =  exp  ~  •  CD3) 


xa  =  fiS.  P3*  . 
a  p+q 


In  the  calculation  of  J2  and  it  is  seen  that  we  are  required  to  evaluate 
the  function  Fm(t),  (nFO,  t  or  2),  defined  as 


1 


F  (t )  =  5  u2|T1  exp  l-tu2)  du 
m  0 


(26) 


Satisfactory  computational  procedures  for  this  auxiliary  function  have  been 
21 

given  by  Shavitt,  and  are  based  on  the  recurrence  relation 


Fm^)  =  tet'F^U)  +  exp  (-t)]/(2m+1)  .  (27) 


It  was  found  that  the  assumption  F^  =  0  was  sufficient  to  assure  eight 

decimal  place  accuracy  for  F0(t),  F,(t)  and  F2(t)  for  0  <  t  £  17.  then  the 

argument  is  larger  than  17  the  asymptotic  forms  given  by  Shavitt  were  employed 

The  remaining  integrals  have  been  treated  within  a  different 

22,23  .  . 

context  and  will  not  be  dealt  with  he'e  in  detail.  There  are,  however, 
several  points  that  should  be  mentioned. 

The  integrals  and  J#  can  be  reduced  to  a  knowledge  of  Dawson's 
Integral,  defined  as 


D(X)  =  exp  (-X2)  5  exp  (y2)  dy 

0 


(28) 


In  order  to  evaluate  this  function  rapidly  and  accurately,  several  different 
representations  were  employed,  depending  upon  the  magnitude  of  the  argu- 

5,26  g  x  <  2j  2  £  X  £  6j  6  <  X).  For  small  values  of  X(<1)  it  wa 


found  convenient  to  use  a  Taylor  series  expansion  of  the  integral  about  X  =  0, 
This  results  in 


D(X)  =  exp  (-X2) 


xan-i 


(2n-1 )  (n-1)J 


(29) 
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The  rapid  convergence  of  the  above  form  is  illustrated  by  noting  that  13  terms 
are  required  for  the  severest  case  (X  =  0.99999999).  Finally,  for  X=0,  which 
represents  a  special  case  in  and  J5,  a  particularly  simple  form  is  found 

J4  =  2ti3/2  Ij/J)1/2  ,  J3  =  2n3*K-yp1/2.q1/2.(p+q)  . 

The  integrals  <4  and  J7  have  been  calculated  numerically.  These  numerical 
integrations  can  be  satisfactorily  computed  by  splitting  the  range  (0,«>) 

27 

into  the  ranges  (0,1)  and  (l,00)  and  using  24-point  Gauss-Lagendre  quadrature 
in  each  range.  A  similar  numerical  technique  was  used  in  the  calculations  of 
'he  electronic  forces.  Here  the  correct  components  [see  equations  (18),  (19) 
and  (20)]  are  obtained  from 

1  63(1) -2-jGc(i)  dv,  (30) 

rA' 

by  th.  identities**1 

exp  U  r^1)  =  ZA,  exp  U  r^)  (31) 

and 

A  exp  (-a  rMJ)  =  VA)  exp  (-a  rA,*)  (32) 

where  now  lA/^3  *n  Equation  (30)  is  replaced  by 


Since  many  of  these  integrals  have  not  been  tabulated  previously,  it 
is  necessary  to  check  the  accuracy  of  the  integral  calculations.  To  accomplish 
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this*  several  identities  can  be  made  that  are  useful.  For  example}  the 
nuR-erical  integration  procedures  that  are  needed  for  <4  were  checked  by 
allowing  nucleus  D  to  become  coincident  with  nucleus  C,  in  which  case  is 
obtained  (this  transformation  does  not  apply  to  HeH+).  A  further  check  on 
the  accuracy  of  the  numerical  integrations  can  be  obtained  by  transforming 
th6  integrals  to  different  forms  before  integrating.  Independent  calculations 
on  the  integrals  -*  J3  can  be  obtained  by  the  use  of  the  Hermitian  property 
of  the  momentum  operator  to  give  c’u'ivalent  forms. 

5.  CONCLUSIONS 

The  rather  large  values  of  the  energy  variance  point  out  the  inability 
of  these  particular  Gaussian  orb.tals  to  describe  accurately  the  true 
wavefunctions  for  the  systems  under  consideration.  While  in  most  cases  the 
electronic  energy  is  predicted  accurately*  the  magnitude  of  U2  is  very  large. 

A  striking  example  is  the  one-electron  H2+  molecule-ion.  In  this  case  the 
energy  is  accurate  to  within  4  x  10"A  hartrees,  while  the  energy  variance  is 
far  from  zero,  .023  hartrees2. 

For  the  two  electron  systems,  any  decrease  in  the  accuracy  of  the 
wavefunctions,  from  the  energy  point  of  view,  is  magnified  in  the  energy 
variance.  This  is  not  unexpected,  and  is  due  primarily  to  the  sensitivity 
of  the  wavefunction  to  the  operators  contained  in  3/2.  Moreover,  the  lack 
of  convergence  of  all  of  the  components  of  2  shows  that  they  each  contribute 
to  these  poor  results.  This  sensitivity  is  even  more  pronounced  for  the 
Gaussian  approximation,  due  to  their  anomalous  behavior  in  those  regions  of 
the  molecule  which  are  heavily  weighted  by  M2. 
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One  obvious  way  to  improve  the  molecular  description  would  be  to  minimize 

the  energy  variance,  rather  than  the  energy,  with  respect  to  the  parameters 

2 

contained  in  the  Gaussian  approximation.  Such  a  calculation  has  the  added 
significance  that,  for  systems  containing  more  than  two  electrons,  not  only 

15 

does  it  require  a  knowledge  of  the  first-  and  second-order  density  matrices 

[see  Equation  (14)  ],  but  also  requires  that  information  concerning  the 

third-order  density  matrix  must  be  obtained.  In  this  way  the  wavef unction 

is  subjected  to  more  rigorous  constraints  which  would  consequently  improve 

its  overall  description.  This  procedure  would  seem  especially  worthwhile, 

since  all  the  integrals  that  arise  can  be  calculated  both  accurately  and 

conveniently.  In  addition,  an  extrapolation  technique,  such  as  the  one 

16 

suggested  by  Conroy,  can  then  be  used  to  give  an  accurate  approximation 
to  the  ground  state  energy.  However,  the  additional  constraints  that  such  a 
calculation  would  impose  on  the  wavef unction  would  necessitate  the  use  of  a 
more  flexible  basis  set.  Investigations  of  this  nature  are  presently  being 


undertaken* 
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C.  G0HNAL  THEORY 

.  INTRODUCTION 

The  need  to  carry  the  accuracy  of  calculation*  on  molecular  *y*ta»s  beyond 
th.t  of  Hartrue-Fock  accuracy1  ha*  led  to  the  development  of  *everal  approache*. 

Among  thee  which  retain  the  physical  Intuition  *uch  a*  that  found  In  Bertree-Fock- 
Roothaan  theory,  the  approach  u.lng  gemlnal  theory  seems  to  hold  some  promt**. 

Although  considerable  atudy  ha.  been  done  on  the  propertle*  of  gemlnal*,  faw 
ab  Initio  calculation*  have  been  done  uslug  the  theory.  One  of  the  reasons 

for  the  paucity  of  *uch  calculations  lie*  In  the  difficulty  in  solving  the  equations 
determining  the  gemlnal*.  Thla  has  led  investigator,  to  the  use  of  simplifying 
approximations  to  produce  manageable  equations18,  or  to  abandon  the  equation,  entirely, 
and  obtain  the  gemiuals  by  an  Infinite  .erie*  of  orthogonal  transformations  of 

the  original  basl*. 

In  this  study,  a  new  procedure  for  the  determination  of  the  beat  separated 
pair  gemiuals  of  a  closed  shell  molecule  is  described.  The  method  treats  the  set 
of  coupled  in te gro-d iff erect iel  equation,  which  determine  the  gemlnal.  without 
simplifying  assumptions,  and  produces  an  algorithm  which  is  quedratically 
convergent  to  the  desired  result. 

2.  GEM  INAL  EC&JATIONS  AND  CONNECTION  TO  HARTREE-fXK-fOOTHAAI'i  THEORY 

The  development  of  the  desired  algorithm  will  be  done, conveniently  by  first 
casting  the  equations  detatminlng  the  gemlnal*  in  a  form  analagous  to  that  found 
in  Hartree-Fock-Roothaan  theory.  To  accomplish  this,  we  recall  that,  for  a  closed 
•hell  syetem  of  gN  electrons,  the  total  electronic  wavefunction  to  be  used  as  e 
trial  wavefunction  can  be  written  in  the  form  of  an  antirysseetrised  product  of 

geainala 

t  (1,£...2N)  .  ft  (Cl  U»2>  =*< 5, 4) *  *  •  Gn( 2M-1, 2M)  J.  (33) 

where  0  is  celled  •  gemlnal,  and  ie  a  function  of  the  space  and  epln  coordinates 
I* 
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of  two  electrons.  The  geminals  are  usually  assumed  to  be  antlsynmetrlc,  l.e., 

G^(l,2)  =  -  0^2,1)  , 

so  that  ^  is  a  partial  antlsymmetrlser  which  exchanges  electrons  between 
geminals  only.  Also,  each  gee  Inal  will  be  written  as  a  product  of  a  space  and 
•pin  part, 

G^l.2)  .  yi,2)  Ml, 2)  ,  (34) 

mnA  in  this  study  we  shall  consider  only  the  case  where  a  singlet  spin  function 


for  each  gsminal  is: 

yi, 2)  .  [  o{l)  0(2)  -  0(1)  or( 2)  )//2  ,  p  «  1,2.. .,N, 

Without  loss  of  generality  ,  the  spatial  part  of  the  geelnals  can  be  written  In 

22 

natural  orbital  fore: 

Au(i, 2)  cpi  Xpl^  V(2)  ,  (35) 

In  which  the  natural  orbitals  are  orthonormal 

I  X*t  (1)  X^U)  <Wi  <=  .  (36) 

and  the  genlnals  are  taken  to  be  weakly  orthogonal,  i.e., 

;jA*(l,2)  Av(l,2)  dVi  dVa  -  .  (37) 

The  strongly  orthogonal  geminals  (or  separated  geminals)  to  be  employed  In  this 


study  have  the  further  resection  that 

f  A*(l,»)  \(l,n)  dVx  =  0  (p*v)  ,  (38) 

whlchk  along  with  Eq.  (5),  Implies  that 

J<i(1)*vJ(1)Wl  'Vv4u  •  (39) 

Using  a  wave  function  of  this  form,  the  energy  expression  can  be  written  as 

cJi  +  V  y  [pi, pJipJ.pl] 

+  £>2  2  c*<  [2(pi,pi|  vj,  vjl  -  [pi,vj|vj,pi)|>  ,  (40) 

v(rfp)  i,j  *  VJ  »■ 
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where 

tnl|h|nO  -Jxj  («  h  JV,  ,  (41) 

-  JJ  X*(l)  Xjd)  lri£  \(2)  Xj(  s)  ,  (42) 

and 

h  »  -  I  »  (43) 

on 


and  the  summation  over  cr  ie  taken  over  a^l  nuclei. 


To  obtain  the  optimum  separated  geminals,  the  energy  expression  of  Equation  (40)  «uat 
be  minimized  simultaneously  with  respect  to  both  the  occupation  numbers  (c  j) 
and  the  natural  orbitals  (x^)*  Independent  variation  of  Equation  (40)  with  raapoct 
to  the  occupation  numbers  yields  the  following  set  of  coupled  eigenvalue  equations 


Hki  clki  "  Vi* 

1=1 


ka  1;2|...  (44) 

ps  1#  2.  •  .,H, 


where 


Hki  “  2C|dt|h||ilel6ki  +  "’^,nl|pl,pk] 


5..?  c*</\[pk,pji|vj,vj]  -fc[pk,vj|vj,pk]\ 

™  W  JJL  v\  J 


(45) 


and  c  Is  the  energy  of  an  electron  pair  in  the  pth  geminal  In  the  average  field  of 

^  IQ 

the  other  electrons.  These  equations  can  be  solved  satisfactorily  hy  choosing 

C  =1  and  =  0  (j>i),  «a  the  initial  guess  in  the  iterative  procedure  to 

obtain  the  solutions  of  Eq.  (44),  and  need  not  concern  us  further. 

Of  greater  interest  are  the  coupled  lntegro-dlfferentlal  equations  that  arise 

from  the  variation  of  Eq.  (40) with  respect  to  the  natural  orbitals.  To  show  the 

connection  to  Hartree-Fock  theory,  as  well  as  to  aid  in  the  later  development  of 


.  r  .. 


Ill 


th«  quadratically  convergent  algorithm,  let  ua  firat  rewrite  the  energy  expression 
of  Eq.  (40) in  matrix  form  aa 


*E  ■  £  *Si  (ill  -  v) +  t  Vi  cuJ  ( iji  v) 

U«1  i,  J 


N 

+  I 

v(*p) 


r  cw  cCi  <2ij  -  £vj>  ^  j 

j 


where  each  natural  orbital  haa  been  written  aa 

v  '  4  V 


•ad  vher.  £  la  a  rov  vector  of  bat  la  orbitals,  and  a  It  a  coliaao  vector  of 


coefficients  needed  to  transform  the  basis  orbitals  Into  the  natural  orbital 

A1*0#  matrices  J^,  and  are  Hermltlan  square  matrices 

with  elements  (M)  ,  given  by 

Pm 


<!!.)„-  ,  UB) 

where  M  Is  a  one  electron  operator  chosen  from  Eq.  (43),  or 

Jvj(1)  *  I  \j(2>  \j(2)  dv*/r**  >  (49) 

or 

KVJ  (X)  X(1)  =  J  ^(1)1  X*j(2)  X  (2)  dVs  #  (50) 

Thus,  the  determination  of  the  optimum  set  of  natural  orbitals  is  accomplished 
by  the  determination  of  the  optimum  vectors  a  .  v 

To  effect  this  determination,  we  find  the  variation  in  the  energy  with  respect  to 
a  variation  in  each  of  the  vectors,  subject  to  the  constraint  that  the  natural 


orbitals  resale  orthonormal.  This  procedure  .exults  In 


6  E  ■  2 


N  Nu  .  /  N  Nv  \ 

rr  •<i  (  vw*r  i  v, J 

1B1  1=1  '  V«1  J»1  ' 

N  N \>  m  f  *  *  N  Nv  *  *  \ 

2  r  r  ( C  -  z  r  s~  vj  w 

U*X  1*1  '  V=1  J*1  • 


where 


F  ,  *  cz. 

1*1  |*1 


a+ 


rvviv j+  £  s 

J»1  v(^i)  j=l 


V  •  «> 


end  S  la  the  overlap  matrix  between  the  basis  orbitals,  with  elements 

<s)ij  -  I  ♦"  <IV  .  (53) 

and  -21^^  ^  la  the  Lagranglan  multiplier  that  haa  been  Introduced  to  assure  the 
orthonoraality  of  x  4  andxvj.  The  condition  that  6E  *  0  for  arbitrary  vectors 
6a^ .  and  6a^.  results  in 

«^Jbi  M^l 


W  * 


*  * 
^1  ^1 


N 

ii 


V*  1 


N  Ny  *  * 

E  E  i  Si* 


wl  J-l 


From  Eqs.  (94)  and  (99)  it  is  clear  that  X  ,  .  Is  an  element  of  a  Hermit tan 

VJ  )|J»1 

matrix  X,  with  X  .  .  «  X  .  ..  Thus,  Eqs.  (54)  and  (55)  s^e  equivalent. 

~  vj ,^i  ,vj 

As  is  apparent  from  the  fora  of  Eq.  ($4),  the  equations  determining  the  best 

natural  orbitals  have  the  same  form  as  the  equations  that  are  encountered  in  multi- 

23 

configuration  Hartree-Fock-Roothan theory.  It  is  not  analagous  to  closed* 

24 

shell  Hartree*Fock*Roothaan  theory,  because  F^  is  not  invariant  with  respect  to 
a  unitary  transformation  of  the  basis  orbitals.  Thus,  the  reduction  of  Eq.  (54) 


113 

to  pseudo-eigenvalue  aquations  by  a  unitary  transformation  of  the  basis  is  not 

applicable  in  this  cast.  Instead,  the  following  quadratically  convergent  algo- 

25 

rlthm  has  been  developed  to  obtain  the  desired  solution  of  Equation  (54) « 


3.  QUADRAT ICALLY  CONVERGENT  PROCEDURE 

Considering  real  orbitals  and  coefficients  henceforth(  it  is  easily  shown 
froa  Iq.  (54)  that 


^pl,vl  ivj  * 

Xvj  ,pi  *  ~^i  ~vj  ~vj 

To  assure  the  preservation  of  synutry,  we  taka 

Xvj,pi  Ji  %i  +  £vj^vj  *  Vi.vJ* 

which  allows  Kq.  (54)  to  be  rewritten  as 

i  2  l4t(v  v  As1 

vml  J*0L 


(56) 

(57) 

(58) 


(59) 


To  establish  the  iterative  procedure,  let  us  assume  that  a  reasonable  approxi¬ 
mation  is  known  for  each  a^  (and  hence  for  each  F^)»  !•••»  and  s^, 
respectively.  Then,  the  desired  solutions  can  be  written  as 


v  ■ 

t 

(60) 

V 

f 

(61) 

and 


where  terns  in  (&?)*  and  (ga^)*  and  higher  powers  have  been  aasuasd  to  be 


11  enough  to  Ignore,  and  the  set  of  vectors  are  chosen  to  oe  orthonormal,  i.e., 

$  ~  iSj  *  V*  •  (62) 

Insertion  of  Iqs.  (60)  and  (61)  into  Kq.  (59)  and  subsequent  neglect  of  tents 
ci  power  two  or  higher  in  (Ag)  or  (AJL)  gives 


F*  a*  .  +  F*.  A*  .  +  AF  a% 
~Vii  *tki  '•'pi  'Ski 


Defining 


Ms  t  £  |  '$ 

vl  J«1  l 


+  *i£i  ^i+  A£vJ^vJ*5vJ 

jf i  - 1/2  S  E  E  ls3  <$» +  Si)’ 

v*l  J*1 


the  use  of  Eqs.  (57),  (64)  and  (65)  allow*  Eq.  (63)  to  be  written  me 

Q£i  “  -Iji^  L£pi  +  ^  V*»vJ  +  A*i»vJ^  A^J  A*1 

v=l  J=1 


■JSi-AJt‘5  2  'S’  n  1 

V=1  v=l  V=1 

(66) 

where  ^  is  formed  by  using  the  approximate  solutions  a?^  e^,  atd 

Fjj  in  Bq.  (53). 

To  express  the  dependence  of  (&F^)  on  the  vector  Aa^,  we  note  that 


H  N 


Vi‘E  Vp*i 


ps=l 

N 

u 

-E 

0p 

p=i 

■*v 

**%i 

■^i+ 

V 

(67) 

(68) 
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Then,  ueing  Bqi.  (48),  (49),  (50),  (52),  and  (68)  we  see  that  the  variation  In 

P 

can  be  written  as 

-  E  vVj  (AM-  Lvj**^"*1*'  <#* 


J^l 
N  N 


E  S’  <1*1,  kttfc  4<^  <!•"»•* + 


v(¥p.)  k=l 

-  -Bfr&ij&lP”*1-']  .  (69) 

Using  this  result,  the  first  part  of  the  last  term  in  Eq.  (66)  can  be  written  as 

h  i  g  'V 

■hi  fii,  vv [v (C.,j ♦ 

N*1  J-l  U 

+  2  2  V  <%°I,vJ+  *.VJ-  iCl.vJ- Cw 

^(*v)  k»l 

•  h  i:  £  &  s  I?  ( vvVv  (&.vj + sSI.vjm 

vwl  J=1  T^l  k*l  ^ 

+  h- VJi*^  (%°i.vi+  2*.^‘  £.w» }  • 

Rearranging  eiwnations  in  Eq,  (70),  and  defining  the  elements  of  a  suparmatrix  M  as 


(70) 


^^i,Hk*  2  ®  X!  ^  |  V^iV[V(^1»vJ+ 

+  C.VJ-  J?l.vi:V  • 


(71) 


Eq.  (70)  can  be  written  as 

hi  i  ‘vVav^j  -  i  f  • 

v*l  J«1  v=l  J«1 

In  a  similar  manner,  the  last  part  of  the  last  term  in  Eq.  (66)  can  be  written  as 


(72) 
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1*2  2  [v  £i]£s  *  f  £ 

V=1  j=d  V=1  J*1 


(73) 


where 


=  2  ~  S  2  £\>J  j  61\vCvj%k^k^i,vJ+  n*-1*^ 

****  **  t  *)  (74) 

+  (l'S)cvJC^[%(2&.vj*  2^i.vj‘  ^i,vj)3/  • 

Finally,  we  note  that  the  other  term  of  intereat  In  Eq.  (66)  can  be  expressed  as 

N  R 


4-V&  ‘  £  £  'Si.vjfcvj  • 

\*=1  J«1 


(75) 


where 


<£°V,vj“  &.VJ*  + 


Collecting  the  results  obtained  thus  far,  we  find  that  Eq.(66)  can  be  written  as 

•  (77) 


where  the  elements  of  the  supermatricee  in  Eq.  (77)  are  given  by 

„o  -O 


(G°)  .  «  6  6 . ,(B°,  -  F°, ) 

V-  jii.vj  p,v  ij  <~p,i 


(78) 


+  £  V*»vJ +  Ai1*^  *  ^i»vj+  i$i.vj+  &.vj 


and 


&  *  it  v  ■  ~  £  z?  ^i.vAj  • 

V=1  J«*l 


(79) 


4.  DISCUSSION 

The  solution  of  Eq.  (77)  for  Qa  provides  an  iterative  procedure  of  the  Newton- 

26 

Raphaon  type,  which  can  be  expected  to  be  quadratically  convergent  to  the  desired 
result  for  a  sufficiently  good  set  of  starting  vectors.  The  analysis  presented 
here  is  similar  in  form  to  a  quadratically  convergent  algorithm  developed  by 
Wesoel^  for  solving  the  general  Hartree-Fock-Roothaan  equations. 

The  present  work  differs  from  Weasel's  in  that  the  equations  arise  from  the  separated 
pair  model,  and  the  dependence  of  the  Fock  matrices  F  ^  upon  the  vectors  £a  .  has 
been  explicitly  evaluated. 
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Considering  the  convergence  properties  of  the  proposed  procedure,  direct 

calculation  of  the  natural  orbitals  of  a  geninal  trial  function  by  this  algorithm 

19-21 

should  provide  a  favorable  alternative  to  the  method  of  2  x  2  rotations  that 

is  often  employed.  Furthermore,  it  has  not  been  necessary  to  neglect  any  off- 

2,18 

diagonal  Lagrangian  multipliers  in  the  above  analysis  in  order  to  obtain 

the  desired  algorithm. 

Concerning  the  use  of  this  algorithm,  several  comments  are  appropriate. 

For  atomic  calculations,  the  coefficients  obtained  from  a  Gram -Schmidt  orthogonal- 

ization  of  the  entire  set  of  orbitals  should  serve  as  suitable  starting  vectors 

o 

a  ..  In  the  case  of  a  molecular  calculation,  the  principal  natural  orbitala  (i.e. , 

-Hi 

the  flrat  natural  orbital  ln  each  geminal  expansion)  may  be  approximated  by 
the  canonical  Hartree-7ock  orbitals,  or  perhaps  by  localized  Hartree-Feck  orbitals. 

It  should  be  noted  that  the  Hartree-Fock  orbitals  may  be  obtained  within  the 
framework  of  this  algorithm  by  simply  terminating  the  natural  orbital  expansion  of 
each  gaminal  to  one  *:erm  and  hence  setting  c  »  1.  In  approximating  the  weakly 
occuplad  natural  orbitals  (x^>  i  >  l)  of  a  molecular  calculation, there  are  several 
possibilities.  One  could  use  the  virtual  orbitals  obtained  from  the  Hartree-Fock 
calculation  as  the  initial  approximation  to  the  weakly  occupied  natural  orbitals. 
Alternatively,  one  could  construct  molecular  orbitals  of  the  proper  symmetry, 
and  orthogonallze  these  orbitals  to  the  principal  natural  orbitals  and  to  each  other. 
Once  the  starting  vectors  are  chosen,  both  the  atomic  and  molecular  calculations 
would  proceed  in  the  same  manner,  with  the  elements  of  G°  and  g°  being  determined. 

A  correction  vector  As  is  determined  from  Eq.  (77),  and  is  added  to  a°,  forming 
a  new  set  of  starting  vectors,  .  These  corrected  vectors  are  then  Gram- 
Schmldt  orthogonalieed  to  insure  their  orthonormality  and  acceptability  at  new 
starting  vectors.  The  iteration  process  is  repeated  until  Aa  -  0. 

Finally,  it  should  be  pointed  out  that  the  procedure  of  solving  Eq<  (77)  will 
give  the  best  set  of  natural  orbitals  for  a  given  set  of  occupation  coefflclants 


c  In  order  to  obtain  the  best  separated  geminal  wavefunetion,  a  macroiteration 

bi¬ 
cycle  should  be  performed,  shuttling  back  and  forth  between  Eqs.  (44)  and  (77). 

One  might  find  it  computationally  convenient  *or  the  G°  matrix  to  be  sym- 

25 

metric.  In  the  work  by  Wessel,  a  symmetric  G°  matrix  was  obtained  by  a  pro¬ 
cedure  which  might  be  employed  here  in  a  similar  manner.  Neglecting  second  order 
terms,  one  has,  from  the  orthonormality  relationship  of  the  vectors. 


a°J  S  Aa 


i  .  +  a0^  S  Aa  ,=0 

pi  ~pi  ~  ~ v j 


from  which  one  may  obtain 


“  T'  ^  1(F°<+  F°4)a°<  *°!  S  Aa  .+  (F°.+  F° .  )a° .  a0*  S  Aa 
2  x  ,  r  .  )  -nii  -A»j  ~  ~ui  ~pi  ~v  j  ~vj  «ui  ~  ** 


v=l  J=1 


which  can  be  rewritten  as 


5ui  ^?ui+  2  2^  ^Ui.vJ  ~\>j  =  °  ’ 


v=l  j=l 

Similar  considerations  concerning  M°.  .,  N°.  . ,  and  D°.  and  addition  to 

-ni.vj  -^i.vj 

Eq.  (78)  of  identities  such  as  Eq.  (82)  are  necessary  if  it  ie  desired  to  com¬ 
plete  the  sytometrization  process  of  G°. 

Applications  of  the  analysis  presented  here  are  presently  underway  on 


various  atomic  and  molecular  systems,  and  will  be  analyzed  and  reported  separately 


in  future  studies. 
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III,  SPECTROSCOP I C  I NVEST IGATIONS  OF  TRANS l ENT  SPEC I ES 

A.  ENERGY  TRANSFER  IN  NAPHTHALENE-B I  ACETYL  MIXTURES 


.  INTRODUCTION 

One  of  the  proposed  methods  of  generating  appreciable  concentrations  of 
excited-state  molecules  is  that  of  photosensitization.  While  our  principal 
aim  is  the  investigation  of  these  excited  species  by  microwave  spectroscopy, 
we  have  found  that  there  is  a  great  need  for  investigation  of  the  primary 
photosensitization  process  and  its  energy  transfer  mechanism  and  efficiency. 
This  last  point  is  particularly  important  for  the  microwave  investigations, 
since  efficient  modes  of  production  of  the  excited  species  are  necessary  in 
order  tc  overcome  detection  sensitivity  problems. 

To  be  more  specific,  the  photosensitization  processes  in  which  we  are 
interested  are  those  in  which  an  electronically  excited  molecule  transfers 
its  energy  to  another  ground  state  molecule  by  1  radiationless  process.  Of 
the  possible  types,  only  those  processes  involving  triplet  state  donors  will 
be  very  effective,  since  singlet  donors  do  not  usually  exist  for  a  sufficient 
length  of  time.  Although  triplet-*singlet  (T-S)  processes  have  been  known 
for  many  years,  the  analogous  triplet+triplet  (T-T)  process  has  been 
established  for  only  a  relatively  short  time.  In  general  the  mechanism  of 
the  transfer  processes  is  poorly  known,  although  several  theoretical 
discussions  have  emphasized  the  importance  of  dipole-dipole  and  exchange 
interactions  between  donor  and  acceptor. 


In  this  work  we  have  investigated  the  gas-phase  transfer  efficiency 
between  triplet  naphthalene  and  singlet  biacetyl  to  form  triplet  biacetyl 
(T-T  process).  The  few  similar  investigations  previously  reported  include 
the  benzene-bi acetyl  system^  and  the  hexaf iuoroacetone-bi acetyl  system.'* 

Some  more  preliminary  studies  of  other  donors  and  acceptors  have  also  been 
undertaken  in  our  laboratory. 

METHOD  OF  ANALYSIS  OF  ENERGY  TRANSFER 

Whilo  several  other  energy  transfer  modes  might  be  considered,  we  limit 

ourselves  here  to  those  shown  in  Figure  II 1-1.  B  indicates  biacetyl,  N 

indicates  naphthalene,  superscripts  give  the  multiplicity  of  the  states,  and 

subscripts  0  and  1  indicate  lowest  and  first  excited  states,  respectively. 

The  straight  lines  are  processes  involving  absorption  or  emission  of 

radiation  while  the  wiggly  lines  are  radiationless  energy  transfer  processes. 

The  3N0  level  lies  at  21,200  cm-1,  the  3B0  at  20,500  cm-1,  and  the  1 B1  at 

6 

about  22,600  cm”1.  The  energy  levels  are  shown  as  single  discrete  states 
although  this  is  not  meant  to  imply  the  absence  of  vibrational  fine 
structure. 

It  has  been  well-established  for  some  time  that  irradiation  of 

0  3 

naphthalene  at  2537  A  leads  to  effective  population  of  the  N0  state. 

Because  of  the  near  resonance  of  the  3B0  and  3N0  levels  it  was  believed  that 

there  should  be  efficient  energy  transfer  from  triplet  naphthalene  to  ground 

state  biacetyl.  This  has  been  observed  by  irradiating  naphthalene-bi acetyl 

o 

mixtures  with  2537  A  light.  Under  these  conditions,  experiments  performed 
with  a  Hilger  medium-quartz  photographic  spectrograph  show  that  there  is 
essentially  no  biacetyi  fluorescence  from  1B1  states,  but  only  phosphorescence 


from  3B0.  Also,  no  naphthalene  phosphorescence  is  observed.  From  these 
results  we  conclude  (1)  that  the  1 Bt  level  is  not  excited,  that  is,  no  S-S 
energy  transfer  occurs,  (2)  that  the  non-radiative  processes  destroying 
3Nq  are  very  effective,  and  (3)  that  3B0  is  produced  from  3NQ  by  T-T  energy 
transfer. 

Based  on  this  mechanism  the  sensitized  phosphorescence  intensity  can 
be  written 


ISP(2537)  =  lNC2537)  •  *|sc  •  (1) 

o 

where  1^(2537)  is  the  light  intensity  absorbed  by  naphthalene  at  2537  A, 

®ISC  *s  ^he  intersystem  ( 1  N-j-* 3NQ )  crossover  efficiency,  and  O^p  is  the 

quantum  yield  of  the  phosphorescence,  that  is,  it  is  the  fraction  of  3N0 

molecules  which  eventually  produce  biacetyl  phosphorescence.  If  now  the  ’B, 

o 

system  is  excited  with  4358  A  radiation,  the  direct  phosphorescence  of 
biacetyl  is  observed.  The  intensity  of  direct  phosphorescence  may  be  written 

Ip (43 58)  =  In(4358)  •  (2) 

where  $p  is  the  fraction  of  1B1  molecules  which  eventually  lead  to 
phosphorescing  3BQ  molecules.  Combining  Equations  (1)  and  (2),  we  get 


/Isp  Y|B(435a)\/  <ip  \ 

\  Ip  AlN(2537>A*|sc>/ 


(3) 


By  Equation  (3) 
value  of  0.15).^ 


it  is  possible  to  calculate  V  relative  to  Cjs  (which  has  a 
The  first  two  ratios  in  Equation  (3)  can  be  measured 
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experimentally)  and  reasonable  estimates  or  limits  can  be  placed  upon  the 
intersystom  crossing  efficiency)  ^  is  eventually  of  interest  to  know 

the  triplet-triplet  transfer  efficiency)  which  can  be  defined  by 
Equation  (4)» 

&SP  =  M 


This  equation  involves  the  assumption  that  all  molecules  which  get  into  the 
1 B,  state  cross  over  into  the  3B0.  Previous  evidence  indicates  that  this  is 
correct  to  within  2$.  Then  Equation  13)  can  be  written 


/isp\  /lB(4358)\  /  1  \ 

=  yp  /  Vn^2^7)/  \4isc / 


15) 


3.  EXPERIMENTAL 

a.  Measurement  fi£  PhcsphOMSCfiUCg  Intensity  Ratios 

Although  it  is  a  task  of  great  difficulty  to  measure  absolute  phosphores¬ 
cence  intensities)  the  measurement  of  relative  phosphorescence  intensities  is 
relatively  easy.  The  principle  components  of  the  instrumentation  used  for 
these  measurements  are  shovn  in  Figure  1 1 1-2.  The  gas  sample  being  irradiated 
is  contained  in  the  pyrex  T-shaped  cell)  C,  which  is  fitted  with  quartz 
windows.  The  long  portion  of  the  cell  is  20  mm  in  diameter  while  the  sidearm 
has  a  diameter  of  10  mm.  This  sidearm  is  positioned  such  that  its  center 
is  Tj-  cm  from  the  entrance  window. 

Irradiation  of  the  samples  is  performed  by  use  of  the  mercury  vapor 
source)  F)  the  monochromator,  E,  and  the  collimator,  D.  Bausch  and  Lomb 
HP-100  and  SP-200  mercury  vapor  sources  have  been  used  for  most  measurements. 


Phosphorescence  apparatus  schematic. 


1.17 


The  monochromator,  E,  is  a  Bausch  and  Lomb  high- intensity  grating  instrument 

fitted  with  a  grating  having  1200  lines/mm.  With  this  grating  the  mono- 

o  o 

chromator  may  be  used  for  irradiation  anywhere  within  the  2000  A  -  7000  A 
region. 

Phosphorescence  (or  fluorescence)  is  observed  at  the  sidearm  of  the 
T-cell  by  the  scanning  monochromator.  A,  a  0.3  meter  McPherson  Model  218 

plane-grating  instrument.  The  1200  line/mm  grating  provides  a  wavelength 

oo  ,  o 

range  of  1050  A  -  10,000  A  with  a  linear  dispersion  of  2o.5  A/mm.  Detection 

is  accomplished  by  the  photomultiplier  assembly,  B,  a  McPherson  Model  650 

unit  fitted  with  an  EMI  9514B  photomultiplier. 

The  high  voltage  for  the  photomultiplier  is  supplied  via  cable  from 

an  Eldorado  Model  201  Uhiversal  Photometer,  which  also  serves  for  measuring 

the  photomultiplier  output  signal  via  cable  a.  The  photometer  output  drives 

a  Bausch  and  Lomb  VOM-6,  2.5  mv  strip-chart  recorder  for  permanent  display 

of  the  emission  signals.  To  decrease  pen  jitter  on  the  sensitive  photometer 

and  recorder  scales,  an  RC  filter  with  a  time  constant  variable  in  the 

range  0-1  second  has  been  attached  at  the  recorder  input. 

In  a  typical  expur iment  on  the  b; acetyl -naphthalene  system  the  biacetyl 

phosphorescence  intensities  are  measured  in  some  convenient  units  (pamps,  say) 

o  o 

with  irradiating  wavelengths  of  2537  A  and  4358  A.  During  these  measurements 
the  geometry  of  the  system  is  unchanged,  and  since  the  photomultiplier  is 
operated  in  its  linear  region  the  phosphorescence  intensity  ratio  Igp/lp 
is  easily  obtained. 

In  most  measurements  that  we  have  performed  the  phosphorescence  intensities 

o 

were  measured  as  the  peak  heights  of  the  3B0**1B0  band  at  "5100  A  as  observed 
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on  the  recorder  tracings.  We  have  also  evaluated  the  intensities  in  some 

cases  by  integrating  over  the  entire  phosphor escence  band.  The  ratios 

I SP/ ' P  evalua*e<*  by  this  method  have  been  found  to  be  within  5-10%  of  those 

obtained  from  peak  heights,  so  the  majority  of  our  measurements  have 

utilized  the  simpler  peak  height  method. 

In  typical  runs  the  phosphorescence  band  is  scanned  several  times  to 

both  higher  and  lower  wavelengths  and  the  average  peak  height  is  obtained 

from  these  scans.  Figure  1 1 1-3  shows  a  typical  pair  of  scans  over  the 
o 

5100  A  biacetyl  phosphorescence  band. 

b.  Measurement  &£  ihfi  Ratio  fi£  Absorbed  Light  Inton&ily 

The  ratio  of  the  number  of  photons  absorbed  at  the  two  wavelengths  has 
been  measured  in  an  indirect  fashion.  What  is  needed  is  1 0(4358)/!^  (2537) 
in  the  region  of  phosphorescence  observed  through  the  sidearm  T  of  the 
sample  cell.  By  straightforward  but  tedious  algebra,  it  can  be  shown  that 

13(4358) 

I N (2537) 

where  A2(B)  =  absorbance  of  biacetyl  at  4358  A  at  some  distance  la,  A^lB)  = 

0 

absorbance  of  biacetyl  at  4358  A  at  some  distance  lj,  A2(N)  and  At(N)  are 

0 

the  same  quantities  for  naphthalene  at  2537  A,  I 0 (4358)  =  the  incident 
0 

intensity  at  4358  A  in  photons/sec,  and  l0(2537)  =  the  incident  intensity 
0 

at  2537  A  in  photons/sec.  For  our  particular  cell  geometry  we  take  1,  =  7  cm 
and  I2  =  8  cm. 


The  absorbances  in  Equation  (6)  have  been  determined  in  the  following 

manner.  A  10  cm  gas  cell  was  filled  in  all  runs  with  a  gas  mixture  that  was 

identical  to  that  used  in  the  T-shaped  phosphorescence  cell.  The  absorbances 

of  the  biacetyl  and  the  naphthalene  were  then  measured  on  a  Bausch  and  lomb 

Spectronic  505  spectrophotometer.  In  these  measurements  it  was  necessary 

o 

to  subtract  the  small  biacetyl  absorbance  at  2537  A  from  the  total  in  order 
to  obtain  the  absorbance  of  naphthalene.^-  Tho  values  obtained  by  this  method 
need  only  be  multiplied  by  7/10  or  8/10  to  give  A,(B  or  N)  and  A2(B  or  N)  of 
Equation  (6). 

Equation  (6)  was  derived  with  the  assumption  that  the  two  components 
strictly  obey  Beer's  Law  individually  and  when  mixed.  Although  some 
deviation  from  Beer's  Law  behavior  is  expected,  our  measurements  and  those 
of  others  have  shown  that  in  the  pressure  range  of  interest  (30  or  40  torr 
total  pressure)  the  deviations  should  produce  errors  of  no  more  than  five 
percent.^ 

The  ratio  of  incident  light  intensities  is  then  the  only  measureable 
quantity  remaining  in  Equation  (6).  This  has  been  measured  by  interposing 
the  detector  element  of  a  Kettering  Model  8301  radiometer  in  the  path  of  the 
light  beam  before  it  enters  the  sample  cell.  A  direct  reading  in 
erg  cm”2  sec-1  is  obtained  for  each  wavelength.  Appropriate  conversion  to 
photons  sec-1  leads  directly  to  the  desired  ratio,  I o (4358/ 1 0 (2537).  Stated 
accuracy  of  the  radiometer  is  5%,  and  the  reproducibility  seems  to  be 
about  2%. 

Some  attempts  were  made  to  calibrate  the  radiometer  by  ferrioxalate 

8 

actinometry.  Since,  however,  the  actinometry  experiments  had  an  uncertainty 
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of  about  10J6,  it  was  possible  to  show  only  that  the  radiometer  was  accurate 
to  within  this  limit. 

In  most  cases  no  attempt  was  made  to  perform  accurate  direct  pressure 
measurements  of  the  biacetyl-naphthalene  mixtures.  Instead,  we  have  used 
the  Spectronic  505  absorbance  measurements  with  the  known  values  of  the 
extinction  coeff icients^*^  to  provide  an  indirect  but  consistent  pressure 
determination  for  each  component. 

4.  RESULTS  FOR  NAPHTHALENE-B I  ACETYL  SYSTEM 

Table  lll-l  gives  the  measured  phosphorescence  ratios  Igp/lp  f°r  seven 
different  runs.  The  value  of  CgCg  is  given  for  each  of  these  runs,  where 
£  g  and  Cg  are  the  extinction  coefficient  and  concentration  in  moles/liter, 
respectively,  of  biacetyl.  With  the  average  value  of  €g  =  10.2,  Cg  and 
hence  the  partial  pressure  of  biacetyl  are  easily  obtained. 

In  Table  lll-l I  we  have  listed  the  experimental  values  of  I (4358)/l (2537) 
for  the  seven  runs  of  Table  lll-l.  The  rather  large  variations  in  these 
values  are  produced  by  two  factors:  (l)  changes  in  incident  light  intensity 
ratios,  and  (2)  changes  in  the  fraction  of  radiation  absorbed  in  the 
viewing  region.  We  found  large  variations  in  the  incident  light  intensity 
ratios  from  day  to  day  and  even  from  hour  to  hour.  However,  in  Rur3  6  and  7, 
the  large  difference  is  caused  by  the  use  of  a  different  SP-200  UV  light 
source.  Also  listed  in  Table  lll-l!  are  (which  is  determined  only 

by  the  vapor  pressure  of  naphthalene  at  room  temperature)  and  1 0<4358)  / 1 0  (2537 )  - 
When  the  experimental  values  of  Tables  lll-l  and  lll-l I  are  substituted 
into  Equation  (3),  along  with  the  value  <pp  =  0.15,  we  obtain  the  values  of 
$Sp$|sc  given  in  Table  Ill-Ill.  The  average  value  has  been  listed  as 
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Table  lll-l.  Ratio  of  Sensitized  to  Direct  Phosphorescence 
Intensities  for  Biacetyl  at  Different  Pressures 


Run 

€bcb 

isp/lp 

1 

.0183 

.0374 

2 

.0079 

.0399 

3 

.0037 

.0762 

4 

.0176 

.0185 

5 

.0147 

.0180 

6 

.0322 

.0263 

7 

.0169 

.0423 

Table 

lll-ll.  Measured  Values  of  the  Ratio  of  the  Number 

of  Photons  Absorbed  at  the  Two  Wavelengths. 

Also 

listed  are€tyCty  and  incident 

light  intensity  ratios. 

Run 

1 o (4358)/ 1 o (2537)  I B (4358)/ 1 N (2537) 

1 

.0123 

3.22 

5.25 

2 

.0160 

8.21 

5.16 

3 

.0214 

8.21 

2.14 

4 

.0153 

9.61 

10.38 

5 

.0143 

11.03 

11.16 

6 

.0142 

0.570 

0.948 

7 

.0125 

0.691 

0.856 
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Table  I 


O 


l-lll.  Sensitized  Phosphorescence  Quantum  Yields 
Times  (Tor  ’N,-*3^). 


%P‘*ISC 


1 

2 

3 

4 

5 

6 
7 


.027 

.030 

.024 

.029 

.030 

.037 

.033 


Average  =  0.030  t  0.003 


mmm 


0,030  ±  0.003,  from  which  it  is  possible  to  calculate  [from  Equation  (4)3 


v4w  -  0.20 


°TT4!3C 


From  the  latter  two  results  we  also  have 


=  6.7 


From  the  results  of  Table  Ill-Ill  we  found  a  statistical  uncertainty 

of  ±0.003.  It  should  be  mentioned  that  several  more  runs  were  performed,  but 

for  any  of  several  reasons  these  results  were  thrown  out.  The  sources  of 

error  which  might  occur  include  (1)  system  geometry  changes,  (2)  impurities, 

(3)  photometer  and/or  radiometer  calibration  errors,  and  (4)  lamp  intensity 

variations.  With  regard  to  (1),  it  was  important  that  no  geometry  changes 

0 

in  the  optical  path  occurred  during  the  irradiations  at  4358  and  2537  A. 

Some  runs  were  thrown  out  for  reasons  attributable  to  this  factor.  The  major 

impurity  that  was  guarded  against  was  oxygen  because  of  its  well-known 

9 

quenching  abilities.  Naphthalene  and  biacetyl  samples  were  well  degassed 
and  were  handled  on  high  vacuum  lines.  In  this  connection  it  is  worth 
mentioning  that  no  appreciable  changes  in  phosphorescence  intensities 
occurred  over  the  period  of  several  hours  so  no  appreciable  photolytic 
decomposition1^  of  biacetyl  was  caused  by  the  rather  mild  irradiations. 

Instrument  calibrations  are  still  of  some  concern  to  us,  but  we  believe 
that  the  cumulative  errors  from  these  factors  could  be  no  more  than  10—1 55t. 

We  mentioned  previously  the  possibility  that  the  peak  height  phosphorescence 
intensities  might  be  in  error.  These  results  were  checked  in  several  runs 
against  the  integrated  intensities.  A  typical  set  of  results  for  two  runs 
are  given  in  Table  1 1 1 —IV. 
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Table  lll-IV.  Integrated  Versus  Peak  Height 
Phosphorescence  Intensity  Ratios 


August  26  Data 

Run  •sp/lp  (Peak)  Igp/lp  (' ntegrated) 


I 

II 

III 


.568 

.185 

.318 


.574 

.201 

.312 


August  27  Data 

I  .374 

II  .290 

III  .257 


.349 

.293 

.243 
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The  problem  which  caused  the  most  difficulty  was  the  variation  of  lamp 
intensity  ratios  with  time.  Many  bad  sets  of  data  were  obtained  and 
discarded  before  the  magnitude  of  this  problem  was  fully  appreciated.  It 
was  found  that  after  the  Bausch  and  Lomb  SP-200  lamp  was  on  for  a  sufficient 
period  of  time  the  intensities  of  2537  and  4358  \  radiation  varied  smoothly 
but  were  not  ever  constant  or  reproducible  from  day  to  day.  Figure  1 1 1-4 
shows  plots  of  these  variations  versus  time  that  were  obtained  on  two 
successive  days  with  the  same  lamp.  Because  of  these  variations  we  took 
great  care  to  perform  the  phosphorescence  measurements  as  rapidly  as  possible. 
In  many  cases  the  lamp  intensity  ratio  was  measured  before  and  after  each 
run  to  insure  that  the  variation  during  the  run  was  not  too  great.  In  any 
case,  we  believe  this  is  probably  the  major  source  of  our  errors,  since  it 
is  difficult  to  be  certain  that  it  has  been  eliminated.  We  hope  to  obtain 
a  more  stable  lamp  so  that  this  problem  may  be  definitely  eliminated. 


5.  DISCUSSION 

In  this  study  it  was  not  possible  to  achieve  concentrations  such  that 
naphthalene  would  be  in  excess,  since  the  naphthalene  pressure  was  limited 
to  the  equilibrium  vapor  pressure  at  ambient  room  temperature  (approximately 
75  microns).  For  this  reason  it  was  not  possible  to  investigate  the  expected 
sigmoidal1  ^  $<,p  versus  pressure  curve.  Hence  our  results  are  best  compared 
to  the  benzene-bi acetyl  results  for  P(biacetyl)/P(C6H6)  »  1.  For  these 
conditions  Ishikawa  and  Noyes  found  ®gp  -0.02.  In  our  work,  ^gp'^igg  =  0*03. 
The  value  of  ®|gQ  is  unknown,  but  has  a  value  of  0.39  in  solution.  If  this 
value  were  valid  in  the  gas  phase  we  would  get  ®^p  =  0*077,  which  is  somewhat 
larger  than  the  benzene-bi acetyl  result  at  comparable  conditions.  We  would 
also  predict  <Pyy  =  0.51  from  our  results. 


Actually,  the  intersyetem  crossing  efficiency  for  benzene  is  larger  in 
the  gas  phase  than  in  solution1,2  so  a  similar  result  would  not  be  unreasonable 
for  naphthalene.  If  this  were  the  case  the  value  of  Gp.  permitted  by  our 
data  would  have  to  decrease,  however.  For  example,  if  =  0.78  (the  value 
for  benzene  in  gas  phase),  our  data  would  yield  ^y  =  0.26  and  =  0.038. 

One  further  means  of  analyzing  our  results  is  to  note  that  for  benzene- 
biacetyl,  -  1.0,  that  is,  every  triplet  benzene  produces  a  triplet 
biacetyl.  If  we  assume  Gy-j-  =  1.0  for  our  system,  we  find  =  0.20  and 
®3P  =  0.15.  Table  lll-V  summarizes  our  experimental  findings  and  the 
resulting  three  possible  sets  of  ®|gQ  and  G^y. 

At  present  our  data  do  not  permit  a  choice  between  the  various  assumptions 
I,  II  and  III  of  Table  lll-V.  Assumption  111  is  probably  the  most  reasonable 
since  at  high  biacetyl/naphthalene  pressure  ratios  the  bimolecular  T-T 

process  should  be  much  more  efficient  than  the  processes  destroying  3No. 

4 

The  results  for  benzene-bi acetyl  are  in  accord  with  this  belief.  Neverthe¬ 
less,  the  results  from  Assumption  II  cannot  be  discounted  even  though  the 
assumption  leads  to  a  rather  small  T-T  transfer  efficiency.  It  is  hoped 
that  further  measurements  will  permit  a  better  distinction  between  the 
various  possibilities. 
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Table  1 


!-V.  Calculated  Values  of  Quantum  Yields 
Based  on  Various  Assumptions. 


_ Experimental 

V^ISC  *TT*ISC 

0.030  0.20 


Calculated 


^ISC 

% 

1 

0.077 

0.39a 

0.51 

II 

0.038 

0.78a 

0.26 

III 

0.15 

0.20 

i.ooJ 

Assumed  value;  see  text. 
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